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ABSTRACT 


By  use  of  computer  simulation,  it  is  possible  to  predict  the  effect 
of  traveling  ionospheric  irregularities  upon  the  data  obtained  by  HF 
sweep-frequency  ground-backscatter  sounders.  The  simulation  process  makes 
use  of  an  ionospheric  model  in  which  radio  rays  are  traced  to  determine 
the  time  distribution  of  energy  along  the  distant  ground  after  skywave 
travel.  Subsequent  computation  yields  a  synthetic  record  of  the  echo 
amplitude  as  a  function  of  radio  frequency  and  time  delay,  much  like  the 
experimental  data.  The  degree  of  realism  achieved  is  limited  primarily 
by  the  need  for  economy;  the  process  used  here  is  a  simplified  version  of 
one  previously  described  (Ref.  1). 

It  is  shown  that  a  localized  ionospheric  irregularity  causes  a  dis¬ 
continuity  in  backscatter  echo  amplitude  with  a  limited  e>tent  in  frequency 
and  with  a  delay  that  is  comparatively  independent  of  frequency.  This  dis¬ 
continuity  has  the  appearance  of  a  short  streak  on  the  synthetic  records; 
such  a  feature  is  common  in  experimental  data.  When  an  ionospheric  irregu¬ 
larity  exists  from  the  ground  to  an  unlimited  altitude  and  is  formed  with 
its  long  axis  tilted  about  45°  away  from  the  observer,  a  backscatter  streak 
appears  which  is  spread  over  the  entire  frequency  range  of  the  background 
echo.  Furthermore,  this  streak  is  tilted  on  the  record,  i.e.,  its  time 
delay  increases  with  an  increase  in  frequency.  When  a  similar  anomaly  is 
programmed  so  that  it  tilts  inward  toward  the  observer,  then  the  resulting 
streak  has  a  delay  that  decreases  with  an  increase  in  frequency.  A  rela¬ 
tion  between  the  tilt  of  the  anomaly  and  the  tilt  of  the  backscatter  streak 
is  thus  demonstrated.  The  relation  is  in  accord  with  a  limited  observation 
that  tilted  streaks  in  experimental  backscatter  move  only  in  the  predicted 
direction  of  anomaly  motion  as  deduced  from  the  measured  direction  of 
anomaly  tilt. 

Anomalies  with  constant  altitude  are  shown  to  produce  new,  short- 
range  leading  edges  on  the  backscatter.  Experimental  data  having  this 
form  can  thus  be  interpreted  as  an  indication  of  the  presence  of  small 
"layerlets"  within  the  ionosphere,  but  it  is  noted  that  there  are  at  least 
three  other  possible  causes  for  this  particular  form  of  backscatter  echo. 
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I  INTRODUCTION 


Many  forms  of  experimental  data  indicate  that  there  are  extensive 
traveling  irregularities  in  the  earth's  upper  atmosphere.  Study  of  the 
mechanics  of  wave  propagation  in  this  regime  indicates  that  a  consider¬ 
able  portion  of  the  traveling  irregularities  are  "atmospheric  gravity 
waves,"  that  is,  waves  whose  propagation  depends  upon  the  force  of  grav¬ 
ity.  In  noctilucent  clouds,  which  are  near  the  lower  edge  of  the  iono¬ 
spheric  K  region,  one  often  sees  a  number  of  long,  straight,  pnrallel 
ripples.  Figure  1  is  an  unusual  set  of  pictures  of  these  clouds.  Here 
we  see  the  intersection  of  two  different  wave  trains  that  are  crossing 
one  another.  Lapse-time  motion  pictures  of  such  clouds  (Ref. 2)  show  that 
the  ripples  move  relative  to  the  apparent  motion  of  the  material  (the 
wind)  and  that  there  may  be  different  ripples  moving  in  different  direc¬ 
tions  simultaneously.  The  existence  of  such  structures  is  also  indicated 
in  the  variable-azimuth  backscntter  data  of  Tvctcn  (Ref.  3),  which  often 
show  the  apparent  progress  of  a  sequence  of  long,  straight  waves  as  they 
travel  through  the  ionosphere.  Similar  irregularities  have  boon  mapped 
in  a  vertical  plane  by  measuring  incoherent  electron  scatter  (Ref.  4)  , 
but  this  experimental  method  works  only  above  the  observer  and  therefore 
is  limited  by  the  immobility  of  the  required  equipment.  Other  experimen¬ 
tal  techniques  can  show  evidence  of  the  passage  of  irregularities,  but 
they  produce  data  from  which  it  is  comparatively  difficult  to  deduce  the 
morphology  of  the  electron  density  distribution. 

It  will  be  shown  here  that  one  can  extract  clues  about  the  vortical 
profile  of  a  distant  traveling  wave  from  sweep-frequency  11 F  sounder  mea¬ 
surements  of  ground  backscntter  which  propagates  each  way  via  the  iono¬ 
sphere.  This  is  not  a  new  form  of  data;  it  is  similar  to  the  data  of 
Hunsuckcr  and  Tvctcn  except  that  here  the  antenna  azimuth  is  fixed  and 
the  frequency  slowly  vnries  across  the  11F  spectrum.  Figure  2  shows  30 
records  of  such  backscntter  obtained  with  a  sounder  at  Stanford.  Notice 
that  a  number  of  near-horizontal  streaks  are  present;  these  enn  be  inter¬ 
preted  to  shed  light  on  the  vertical  distribution  of  electrons  within 
traveling  ionospheric  irregularities.  Since  skywnve  energy  travels  at 
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a.  Photograph  taken  by  K.  A.  Devine  at  0115  AST 
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b.  Photograph  taken  by  K.  A.  Devine  at  0133  AST 


Fig.  1.  N0CT1 LUCENT  CLOUDS  PHOTOGRAPHED 
SEL-68-029  2 


FROM  FORT  CHIMO  ON  JULY  6,  196-1. 
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c.  Photograph  taken  by  J.  L.  Veillcux  at  0130  AST 
Fig.  1.  CONTINUKD. 


nearly  the  speed  of  light,  the  ground  from  which  this  scatter  is  obtained 
is  approximately  1500  km  distant  for  each  10  msec  of  group  delay.  It 
will  be  shown  hero  that  an  ionospheric  irregularity  tends  to  focus  energy 
at  the  ground  range  of  those  rays  that  have  their  apogee  located  in  or 
near  the  irregularity.  Consequently,  for  one-hop  backscatter,  the  range 
of  a  detected  irregularity  is  approximately  half  the  range  of  the  scat¬ 
tering  ground,  or  750  km  per  10  msec  of  backscatter  group  delay.  By  this 
approximate  rule,  it  may  be  seen  that  the  irregularities  shown  on  the  one- 
hop  backscatter  of  Fig.  2  arc  locatoa  in  the  ionosphere  at  a  distance  of 
600  to  1800  km  from  the  sounder.  (Similar  streaks  in  multihop  backseat- 
tor  will  not  be  discussed  horo,  although  they  arc  plentiful  on  the  data.) 
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Data  of  7  Nov  64 
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Fig.  2.  THIRTY  EXPERIMENTAL  SWEEP-FREQUENCY  BACKSCATTER  SOUNDINGS 
OBTAINED  WITH  A  5  kW  SOUNDFR  NEAR  STANFORD.  (Ref.  5) 
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Fig.  2.  CONTINUED 
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On  September  19,  1967  my  colleague,  Mr.  Taylor  Washburn,  ran  a  tost 
to  check  for  the  existence  of  moving  streaks  on  sweep- frequency  backscat- 
ter.  He  obtained  58  frames  of  backscatter  in  a  time  so  short  that  it  is 
possible  to  observe  the  motlor  of  individual  streaks.  For  ease  of  refer¬ 
ence,  these  frames  are  considered  to  constiti  to  Fig.  3,  and  they  are  pre¬ 
sented  as  an  appendix  at  the  end  of  this  report.  It  is  intended  that  the 
re  ider  should  be  able  to  flip  through  the  appendix  with  his  gaze  fixed  on 
each  frame  to  see  the  motion  of  the  backscatter  irregularities. 

On  the  right  side  of  each  frame  is  given  the  Pacific  daylight  time. 
This  is  Greenwich  mean  time  minus  7  hours.  There  was  nominally  one  frame 
every  2  min  obtained  with  a  sweep  rate  of  250  kHz  per  second  downward  from 
27  to  9  MHz.  The  chirp  sweep  technique  was  used  with  a  transmitted  power 
of  30  kW  into  an  antenna  recently  built  by  IT'  at  the  Lost  Hills  field 
site  of  Stanford  University.  This  antenna  is  a  broadside  array  of  folded 
tilted  monopoles  directed  along  a  bearing  of  86°.  The  antenna  has  18  ele¬ 
ments  on  nominal  40-ft  spacings  for  a  total  aperture  of  680  ft  when  it  is 
in  the  transmission  mode.  The  backscatter  was  received  at  Site  514  on  the 
Stanford  campus  by  means  of  a  rhombic  antenna.  These  particular  data  were 
obtained  during  an  early  test  of  some  new  equipment  and  therefore  many  of 
the  settings  of  equipment  variables  were  not  optimized.  As  a  consequence, 
it  appears  that  the  data  have  a  very  limited  dynamic  range,  which  is  par¬ 
ticularly  unfortunate  in  the  present  context  since  interest  here  is  cen¬ 
tered  on  very  slight  discontinuities  in  the  echo  amplitude.  Many  of  the 
tilted  streaks  on  the  backscatter  are  only  barely  visible  for  a  single 
frame;  it  is  therefore  not  possible  to  follow  the  motion  of  very  many  of 
the  interesting  features.  Doubtlessly,  we  will  obtain  much  better  data 
of  this  type  in  the  near  future,  but  at  the  time  of  this  writing  (Decem¬ 
ber,  1967)  these  data  were  the  best  available  of  this  type. 
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II  THE  METHOD  OF  SIMULATING  BACKSCATTER 


Analysis  of  the  phenomena  that  play  a  part  in  the  production  of 
ground  backscatter  shows  that  each  mechanism  is  comparatively  simple 
but  that  there  is  an  interplay  of  such  a  variety  of  mechanisms  that  the 
entire  process  is  rather  complex.  Fortunately,  since  computers  share 
this  characteristic,  backscatter  can  be  simulated  in  a  straightforward 
manner  :  It  is  only  necessary  to  simulate  each  mechanism  and  to  program 
an  appropriate  interaction  of  the  parts,  copying  the  natural  process. 
Because  of  economic  considerations,  it  is  convenient  to  break  the  calcu¬ 
lations  into  two  parts,  designated  "raytracing"  and  "backscatter  synthe¬ 
sis,"  which  are  related  to  one  another,  as  shown  in  Fig.  4. 


Ionospheric 

Model 

♦ 


Ray  Trace 
Program 


/ 

Ray 

Drawings 


Ray  sets 
(Cards) 


Other 

Experimental 

Conditions 


♦ 

Synthetic  Backscatter 
Plats 


Fig.  4.  THE  SIMULATION  PROCEDURE. 


A.  Raytracing  and  the  Rayset  Method 

To  begin  the  process,  one  must  select  and  encode  an  ionospheric 
model,  which  is  a  two-dimensional  electron  distribution.  Radio  rays  are 
then  traced  in  the  model  by  any  one  of  a  variety  of  methods.  These  will 
not  be  discussed  here  because  they  constitute  an  extensive  subject  in 
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themselves  and  because  the  specific  choice  oi  a  method  is  unimportant  in 

-K- 

the  present  context.  The  program  used  here  is  a  repeated  application  of 
Snell's  law  in  a  refractive  index  distribution  which  is  scalar  because 
the  geomagnetic  field  is  neglected.  The  computed  output  consists  of  ray 
drawings  and  also  of  ray  information  coded  into  punched  cards,  which  are 
called  raysets.  These  contain  enough  information  to  describe  the  time 
and  space  distribution  of  an  ionospherically  propagated  (skywave)  pulse 
along  the  ground.  Thus  the  raysets  provide  the  backscatter  synthesis 
program  with  all  necessary  information  relating  to  the  effect  of  the 
ionospheric  model. 

As  a  consequence  of  this  approach,  it  is  possible  to  compute  the 
effect  on  backscatter  of  any  ionospheric  irregularity  through  which  one 
can  trace  rays.  This  is  the  key  to  the  power  of  the  technique.  With 
the  present  state  of  the  art,  our  ability  to  trace  rays  is  limited  pri¬ 
marily  by  the  necessity  for  economy;  the  methods  fo^  simulation  are  well 
developed. 

B.  Backscatter  Synthesis 

The  synthesis  can  be  carried  out  with  striking  realism;  it  can  be 
compared  to  complicated  bookkeeping  in  the  sense  that  each  simulated 
mechanism  is  simple  but  many  of  them  interact  so  that  the  entire  process 
is  complex.  The  most  realistic  simulation  (Refs.  1  and  6)  required  rough¬ 
ly  1  min  on  a  large  computer.  For  this  study  it  was  desired  to  compute 
backscatter  for  10  to  30  frequencies  per  "sweep-frequency"  record,  but 
the  price  of  10  to  30  min  of  computer  time  was  prohibitive.  The  original 
simulation  program  made  use  of  a  number  of  approximations:  For  example, 
the  geomagnetic  field  was  neglected  during  ray  calculation  and  the  iono¬ 
spheric  model  was  only  two-dimensional.  (Backscatter  has  not  yet  been 
simulated  in  this  manner  without  those  two  assumptions.)  To  reduce  costs 
further,  a  number  of  additional  simplifying  assumptions  have  been  made 
here,  as  follows: 

* 

The  interested  reader  can  refer,  for  example,  to  the  January  1968  issue 
of  Radio  Science ,  which  is  unvoted  exclusively  to  raytracing.  Also,  the 
author  is  preparing  Technical  Report  No.  112,  which  will  describe  the 
Stanford  University  activities  in  high  frequency  radio  raytracing. 
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(1)  It  is  assumed  that  rays  that  travel  out  to  the  scatterer  via 
a  particular  path  can  return  only  by  the  same  path. 

(2)  Only  one  ionospheric  hop  is  allowed  along  the  path. 

(3)  The  ground  scatter  coefficient  varies  with  incidence  angle  and 
radio  frequency,  but  it  no  longer  is  a  function  of  range. 

(4)  For  the  purpose  of  calculating  absorption,  a  single  value  of 
solar  zenith  angle  is  used,  independent  of  range. 

(5)  The  effect  of  pulse  duration  and  shape  is  not  calculated. 

(6)  Constant  factors  are  ignored  since  only  relative  echo  amplitude 
is  calculated. 

With  these  six  additional  assumptions,  the  simulation  program  became 
simplified  to  such  an  extent  that  program  execution  was  reduced  to  rough¬ 
ly  2  sec  per  frequency  sample  in  the  synthetic  backscatter.  This  makes 
the  newer  method  practical  as  a  research  tool  for  preliminary  studies  and 
as  a  means  for  seeking  gross  behavior  patterns.  The  original  method  can 
still  be  used  for  the  study  of  effects  that  would  be  absent  from  this  fas¬ 
ter  program  due  to  the  loss  of  realism.  In  both  programs,  the  sounder 
antenna  has  gain  which  is  a  function  of  both  elevation  and  azimuth  and 
which  furthermore  is  different  at  each  radio  frequency.  It  was  believed 
that  the  gain  pattern  has  such  a  strong  influence  on  the  backscatter  that 
its  form  could  not  be  neglected. 

The  new  approximations,  applied  to  Eq.  (8)  from  Ref.  1,  lead  to  the 
following  expression  for  the  backscatter  power  conveyed  to  and  from  the 
scattering  ground  between  a  neighboring  pair  of  closely  spaced,  precom¬ 
puted  rays: 


Power  = 


G  <3)  aQ(^) 


/ cos  g  \2 
\sin  ' 


(A3)2 

AT  AR 


where 


G  =  antenna  gain,  a  function  of  takeoff  angle 
0  =  ray  takeoff  angle  measured  relative  to  the  horizon 
^  =  ray  landing  angle  at  the  scatterer  end  of  the  path 
T  =  group  delay  along  the  ray 
R  =  ground  range  traversed  by  a  ray 


A3 1  AT,  AR  =  differences  of  the  parameters  for  two  adjacent  rays 
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\|r,  R  =  averages  of  the  parameters  for  the  two  rays 

nQ  =  ground  scatter  coefficient,  set  equal  to  sin*5  ^  with 
optional  constant  b 

The  power  is  reduced  by  absorption,  which  is  calculated  from  the 
formula: 


1700 (cos  0.881  X)1’3 


(F  +  1.2  MHz)' 


Vl  -  0. 


dB  round  trip 


(2) 


2  - 
97  cos  (3 


This  formula  has  been  derived  from  Laitinen  and  Haydon  (Ref.  7)  with 
provisions  that  the  sunspot  number  is  100  and  the  electron  gyrof requency 
is  1.2  MHz.  Here  X  is  the  average  solar  zenith  angle  along  the  path, 
a  single  constant  set  in  the  program  by  the  user. 

The  baekscatter  simulation  process  is  now  very  simple.  At  each 
frequency,  each  neighboring  pair  of  rays  is  used  to  calculate  the  power 
that  would  be  received  according  to  Eqs.  (1)  and  (2).  The  power  is  as¬ 
sumed  to  be  received  uniformly  during  the  time  interval  from  the  minimum 
round  trip  time  to  the  maximum  round  trip  time  (the  minimum  an  maximum 
of  2T  from  the  two  rays).  At  any  moment  in  time,  one  may  receive  the 
sum  of  several  different  power  increments  from  different  radio  propaga¬ 
tion  modes.  The  computer  adds  all  possible  power  increments  as  a  func¬ 
tion  of  time  at  each  frequency  and  generates  a  plot  of  power  vs  time  at 
that  frequency. 

C.  The  Data  Display 

By  trial,  it  has  been  found  most  convenient  to  plot  the  logarithm 
of  power,  normalized  so  that  a  specified  total  number  of  decibels  of  vari¬ 
ation  is  plotted.  This  specification  is  called  "dynamic  range"  because 
it  controls  the  appearance  of  the  data  much  as  the  dynan ic  range  of  a 
receiver  controls  the  appearance  of  actual  baekscatter.  The  data  pre¬ 
sented  here  have  all  been  calculated  at  30  dB  to  provide  uniformity, 
although  the  decibel  rpecif ication  is  arbitrary. 

The  display  of  the  synthetic  data  is  designed  to  simulate  the  ex¬ 
perimental  display  such  as  that  used  on  Fig.  2.  A  difficulty  arises 
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because  the  data  are  three-dimensional.  The  experimental  records  show 
amplitude  by  varying  the  oscilloscope  light  intensity  to  be  photographed, 
but  no  computer  plotting  device  available  at  Stanford  could  produce  a 
variable  intensity  plot.  To  overcome  this  limitation,  the  abscissa  of 
the  synthetic  display  is  used  for  two  different  parameters  according  to 
the  logic  illustrated  on  Fig.  5.  At  the  top  is  shown  a  model  with  three 
dimensions:  time,  frequency,  and  power.  The  synthesis  is  essentially  a 

fixed-frequency  process  that  is  applied  repeatedly  at  multiple  frequencies. 
(The  same  is  true  of  the  experimental  process  when  a  pulsed  sounder  is 
used.)  If  the  fixed-frequency  records  of  Fig.  5a  were  hinged  about  their 
baseline  so  that  they  could  be  folded  down  onto  the  frequency-time  plane, 
the  record  would  look  like  that  of  Fig.  5b.  It  is  necessary  to  adjust  the 
power  variation  scale  so  that  no  single-frequency  record  overlaps  the  base¬ 
line  of  the  next  highc  -  frequency.  With  this  precaution,  the  abscissa  can 
be  used  to  represent  both  frequency  and  power  in  an  unambiguous  manner. 
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g.  5.  THREE-DIMENSIONAL  DATA.  (a)  Presented  in  a  two-dimensional  display; 
(b;  where  the  abscissa  represents  both  frequency  and  power. 


Ill  FIVE  MODELS  OF  THE  IONOSPHERE 


Here,  five  ionospheric  models  will  be  used  to  produce  synthetic 
sweep-frequency  backscatter  records.  One  model  is  the  undisturbed  am¬ 
bient,  a  single  Chapman  layer  with  a  maximum  density  of  10^  electrons/ 

3 

cm  at  an  altitude  of  300  km  and  a  scale  height  of  100  km.  The  remain¬ 
ing  models  are  identical  except  that  each  has  an  anomaly  of  some  type. 

Three  models  contain  a  localized  irregularity  approximately  1200  km 
from  the  hypothetical  radar,  with  shapes  as  illustrated  on  Fig.  6.  The 
top  model  (a)  will  be  called  the  "local"  model  since  it  is  concentrated 
in  a  circular  area.  The  second  model  (b)  will  be  called  an  "outward  tilt" 
because  the  structure  tilts  away  from  the  sounder.  The  third  model  (c) 
will  he  called  an  "inward  tilt"  since  it  is  tilted  inward.  In  each  of 
these  three  models,  the  electron  density  profile  across  the  local  dis¬ 
turbance  has  a  depletion  region  in  which  the  number  density  vs  distance 
is  parabolic  and  in  which  the  central  density  is  30  percent  less  than  the 
ambient  density  outside  the  irregularity  at  that  height. 

The  decrease  of  electron  density  by  30  percent  of  ambient  in  a  trav¬ 
eling  disturbance  is  not  unreasonable.  Reference  8  shows  records  of  a 
disturbance  in  which  the  overhead  electron  density  decreased  by  more  than 
70  percent  within  a  few  minutes  during  the  passage  of  an  irregularity. 
Also,  one  of  the  two  profiles  given  by  Ref.  4  shows  a  large  region  in 
which  electron  density  is  decreased  by  30  percent. 

The  fifth  model  contains  a  small  concentric  "layerlet,"  that  is,  a 
region  of  limited  depth  within  which  the  electron  density  is  perturbed 
in  a  manner  that  is  a  function  only  of  altitude.  For  this  layer,  the 
entire  electron  density  model  can  be  described  by  the  single  curve  of 
electron  density  vs  height  shown  in  Fig.  7.  The  model  was  derived  by 
multiplying  the  ambient  Chapman  function  by  a  single  cycle  sinusoid  of 
density  vs  altitude.  The  maximum  density  change  was  10  percent  and  the 
cycle  wavelength  was  20  km  in  depth,  as  shown  in  Fig.  6d.  In  effect,  we 
have  moved  roughly  6  percent  of  the  electrons  from  the  250-  to  260-km 
interval  up  to  the  260-  to  270-km  interval.  (This  latter  view  is  an  ap¬ 
proximation  since  the  total  number  of  electrons  actually  increased  slight¬ 
ly  due  to  multiplication  b\  the  sinusoid.) 
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It  is  somewhat  debatable  that  this  model  should  be  Included  here  In 
a  discussion  of  traveling  irregularities.  The  multiplying  function  might 
represent  a  traveling  wave  that  is  moving  up  or  down  while  remaining  con¬ 
centric  with  the  earth,  but  it  might  also  represent  a  stationary  shelf  of 
ionization  somewhat  like  a  miniature  ionospheric  luyer.  The  experimental 
data  available  indicate  the  existence  of  these  structures  but  do  not 
provide  clues  about  possible  motion. 
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Fig.  6.  VERTICAL  PROFILES  OF  THE  FOUR  IONOSPHERIC  ANOMALY  MODELS. 
These  are  used  as  perturbations  in  a  single  Chapman  layer  model. 
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Fig.  7.  THE  CONCENTRIC  ANOMALY  USED  FOR  SIMULATION.  (Same  as  Fig.  6d) 


IV  SYNTHETIC  BACKSCATTER  FOR  THE  FIVE  MODELS 


A.  Backscatter  with  the  Ambient  Ionosphere 

The  ambient  ionospheric  model  was  used  to  produce  the  synthetic 
backscatter  shown  on  Fig.  5.  Since  the  model  is  smooth,  it  follows 
that  all  of  the  small  detail  on  the  automatically  drawn  figure  is  at¬ 
tributable  to  errors  in  the  simulation  process.  Some  small-scale  jit¬ 
ter  is  introduced  b>  the  automatic  plottei ,  but  most  is  due  to  the 
introduction  of  numerical  errors  through  the  use  of  finite  steps  in 
the  simulation  process.  In  a  way,  this  is  the  synthetic  analogy  of 
noise  and  it  can  be  reduced  through  the  use  of  smaller  stops  in  tho 
computer  with  a  consequent  increase  in  the  cost.  One  must  therefore 
compromise  and  choose  computation  steps  largo  enough  to  be  economical 
and  yet  small  enough  that  the  desired  results  may  be  percoived.  This 
is  much  like  tho  choice  of  signal-to-noise  ratio:  it  would  be  extrav¬ 
agant  to  strive  for  too  much,  but  the  ratio  must  be  high  enough  that 
the  signal  can  bo  identified. 

Since  each  of  the  remaining  ionospheric  models  is  the  same  ambient 
Chapman  layer  with  a  single  added  anomaly,  we  may  compare  the  resulting 
synthetic  backscatter  to  that  shown  on  Fig.  5  and  attribute  any  differ¬ 
ences  to  the  anomalies.  In  this  way,  the  effect  of  each  anomaly  can  bo 
separated  from  that  duo  to  other  mechanisms  which  influence  the  time  and 
frequency  distribution  of  backscatter  power. 

B.  Backscatter  with  Local  or  Tilted  Anomalies 


The  first  variation  to  be  studied  is  the  effect  of  the  "local"  anom¬ 
aly  in  which  minimum  electron  density  occurs  at  tho  centerline  1200  km 
from  the  simulated  HF  sounder.  The  electron  density  as  a  function  of 
height  at  the  centerline  is  a  parabola  with  its  minimum  at  200  km.  Tho 
outer  boundary  of  the  depletion  region  would  be  a  100-km  circle  if  plot¬ 
ted  in  cartesian  coordinates  of  height  vs  ground  range,  but  because  of 
geocentric  curvature,  tho  actual  boundary  departs  slightly  from  a  circle. 

For  comparison  with  subsequent  models,  the  most  important  feature 
of  tho  local  anomaly  is  its  limitation  in  height:  There  is  no  change 
from  ambient  below  150  km  or  above  250  km.  Tho  effect  of  this  spherical 
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depletion  on  radio  rays  is  similar  to  tho  effect  of  a  magnifying  lens  on 
light  rays;  that  is,  the  rays  that  are  off-center  are  refracted  toward 
tho  conterline.  The  main  complication  arises  here  because  the  conver¬ 
gence  effect  is  mixed  with  ionospheric  refraction,  making  it  somewhat 
difficult  to  visualize  the  total  effect. 

Synthetic  backscattor  computed  with  the  local  anomaly  is  shown  in 
Fig.  8a.  When  compared  to  that  dorived  from  the  ambient  model,  it  is 
seen  that  there  is  a  new  "spike"  in  the  backscattor  at  a  time  delay  of 
about  16  msec,  which  is  the  round  trip  delay  to  the  ground  at  twice  tho 
range  of  the  anomaly  itself.  Ray  plots  show  that  there  is  focusing  at 
2400  km  on  the  ground,  as  expected.  It  is  particularly  significant  that 
this  spike  exists  only  over  a  limited  range  of  frequencies  because  of 
the  limitoH  altitude  interval  occupied  by  the  anomaly.  At  low  frequen¬ 
cies,  none  of  the  rays  is  able  to  penetrate  the  ionosphere  to  the  level 
of  tho  bottom  of  tho  anomaly.  When  the  frequency  is  too  high,  all  the 
rays  go  over  the  anomaly.  Thus  it  is  seen  that  the  altitude  characteris¬ 
tic  of  the  anomaly  influences  the  frequency  characteristic  of  the  back- 
scatter  data.  We  will  see  that  this  altitude-frequency  correspondence 
can  be  carried  a  step  further  to  a  more  useful  cause-effect  relationship. 

The  "outward  tilt"  is  the  next  anomaly  inserted  in  the  Chapman  iono¬ 
sphere  and  its  effect  is  shown  in  the  backscatter  of  Fig.  8b.  This  de¬ 
pletion  region  extends  from  the  ground  upward  and,  because  of  its  unlimited 
altitude ,  tho  backscattor  is  influenced  over  a  wider  range  of  frequencies. 
Perhaps  the  most  important  point  to  notice  is  tho  increase  in  the  delay  of 
the  backscattor  spike  with  an  increase  in  frequency.  This  is  the  key  fea¬ 
ture,  indicating  that  the  anomaly  was  tilted  outward  relative  to  the 
sounder. 

Tho  use  of  the  inward  tilt  produces  the  backscattor  of  Fig.  8c,  again 
unlimited  in  frequency  because  the  anomaly  is  unlimited  in  height.  Now, 
however,  the  delay  of  the  backscatter  spike  decreases  with  increasing  fre¬ 
quency,  indicating  that  the  anomaly  is  tilted  inward.  The  magnitudes  of 
the  tilt  angles  such  as  those  given  on  Fig.  8  are  of  course  dependent  on 
the  choice  of  abscissa  and  ordinate  scales,  so  the  t.'.lt  changes  would  be 
more  apparent  if  the  display  were  compressed  horizontally.  In  a  previous 
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Fig.  8.  SYNTHETIC  BACKSCATTER  COMPUTED  WITH  THE  VARIOUS  ANOMALIES. 
(Those  should  be  compared  to  Fig.  5.) 


21 


SFL-6 8-029 


study,  the  author  used  a  vertical  anomaly  that  produced  ray  focusing  on 
the  ground  at  a  fixed  range  independent  of  frequency  (unpublished).  That 
result  motivated  the  present  work,  since  it  led  to  the  suspicion  that  a 
tilted  anomaly  would  produce  a  backscatter  striation  whose  delay  was  a 
function  oi  frequency. 

The  fact  that  the  tilt  of  an  anomaly  affects  the  tilt  of  the  asso¬ 
ciated  streak  is  perhaps  the  most  important  conclusion  of  this  study. 

On  experimental  backscatter  data,  tilted  streaks  or  striations  are  often 
seen.  One  set  of  data  is  available  with  repetitive  backscatter  frames 
taken  for  2  hours  by  a  fast  sounder  so  that  it  was  possible  t:>  follow  the 
progress  of  individual  striations.  This  is  given  as  Fig.  3  (Appendix). 

Six  moving,  tilted  streaks  can  be  found  on  the  figure.  In  each  of  these 
examples,  the  delay  decreases  with  increasing  frequency  and  the  entire 
echo  moves  inward  toward  the  sounder.  This  echo  behavior  corresponds  to 
an  inward  tilt  moving  toward  the  sounder;  the  direction  of  travel  thus 
agrees  with  theoretical  predictions  that,  in  the  F  region,  tiited  gravity 
waves  usually  move  only  in  the  direction  of  tilt  (Ref.  9).  While  these 
examples  are  not  sufficient  to  support  unambiguous  conclusions,  it  is 
nevertheless  encouraging  that  the  direction  of  motion  and  the  tilt  direc¬ 
tion  agree  with  theory  in  each  case.  Further  experimental  investigation 
is  clearly  desirable.  Also,  we  must  investigate  to  see  if  other  iono¬ 
spheric  models  can  cause  tilted  striations,  for  there  is  no  guarantee 
that  the  interpretation  given  here  is  unique. 

C.  Backscatter  with  a  Concentric  Anomaly 

When  the  concentric  anomaly  (actually  a  thin  ionospheric  layer)  is 
inserted  between  250  and  270  km,  the  calculated  backscatter  acquires  the 
appearance  shown  in  Fig.  8d.  Comparison  of  this  record  with  the  ambient 
standard  shows  a  new  leading  edge  of  the  ground  scatter  with  a  shortened 
delay.  It  might  be  said  that  there  is  now  a  false  leading  edge — false  in 
the  sense  that  it  represents  the  effect  of  a  small  perturbation  on  the 
larger  ionospheric  layer  which  alone  would  produce  only  the  second  (higher) 
leading  edge  of  the  two  shown  on  the  figure. 

It  was  both  gratifying  and  disturbing  to  find  this  effect.  The  dis¬ 
covery  v as  gratifying  because  it  permits  the  explanation  of  multiple,  weak 
leading  edges  that  are  often  seen  in  sweep-frequency  backscatter.  Exam¬ 
ples  are  shown  in  Fig.  2,  frames  (d) ,  (e) ,  (g),  (h) ,  (j) ,  (k) ,  (m) ,  (q) , 
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(s)  ,  (t;,  (u)  ,  (v)  ,  (a),  (3),  (&)  ,  (p,  (it),  and  (n).  Frame  (m)  is  prob¬ 
ably  the  best  example.  It  thus  becomes  possible  to  explain  these  common 
features  by  attributing  them  to  small  ledges  in  the  ionization  profile. 

One  would  expect  that  such  ledges  should  be  common  since  measures  of  elec¬ 
tron  density  vs  height  acquired  by  radio  sounding  or  by  rocket  probes  usu¬ 
ally  show  a  profile  much  rougher  than  the  idealized  Chapman  layer. 

The  phenomenon  is  somewhat  disturbing  if  one  anticipates  using  back- 
scatter  sounding  as  a  means  for  determining  the  ionospheric  structure  at 
a  distance.  The  difficulty  arises  because  multiple  leading  edges  can  be 
caused  by  this  mechanism  and  also  by  a  number  of  others.  Consequently, 
the  appearance  of  a  multiple  leading  edge  may  not  be  a  unique  indicator 
of  a  particular  ionospheric  structural  form,  but  it  is  hoped  that  analy¬ 
sis  will  show  how  to  discriminate  among  the  possible  causes.  Multiple 
leading  edges  can  be  caused  by  magnetoionic  splitting  since  the  ordinary 
and  extraordinary  ray  families  have  different  skip  distances.  They  may 
also  be  caused  by  the  presence  of  different  ionospheric  structures  on  dif¬ 
ferent  azimuths.  These  differences  affect  the  sounding  due  to  the  almost 
unavoidable  presence  oi  sidelobes  and  backlobes  in  the  antenna  beam. 

Further,  the  double  leading  edge  can  be  caused  by  the  separate  action 
of  time  focusing  and  range  focusing,  a  factor  that  has  apparently  not  been 
recognized  before.  To  see  what  is  meant  here,  refer  to  Eq.  (1)  and  notice 
that  the  denominator  contains  both  AT  and  AR.  This  means  that  the  power 
received  will  be  very  large  when  either  the  time  delay  between  adjacent 
rays  approaches  zero  or  the  range  difference  between  adjacent  rays  ap¬ 
proaches  zero.  It  has  long  been  known  that  AR  approaches  zero  at  the 
skip  distance,  whereas  AT  approaches  zero  at  the  minimum  time  delay  and 
the  two  do  not  occur  simultaneously  (Ref.  10).  At  any  fixed  frequency, 
the  first  backscatter  to  be  received  will  be  that  from  minimum  time;  Peter¬ 
son  has  pointed  out  that  this  onset  of  the  signal  is  comparatively  strong 
because  of  the  tine  focusing,  here  indicated  by  AT  in  the  denominator  of 
Eq.  (1).  In  addition,  when  the  energy  returns  later  from  the  skip  distance 
where  AR  approaches  zero,  there  should  be  another  signil  enhancement  at¬ 
tributable  to  the  focusing  of  the  rays  in  space.  It  might  be  said  that  the 
combined  action  of  the  antenna  and  the  ionosphere  produces  an  extremely 
high  gain  at  this  particular  distance. 

Because  there  are  two  separate  sources  of  focusing,  spatial  and  tem¬ 
poral,  it  would  be  expected  that  the  leading  edge  would  be  double.  This 
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effect,  has  not  been  experimentally  observed,  partly  because  the  time  sep¬ 
aration  between  the  enhancements  is  small,  partly  because  the  effect  is 
mixed  in  with  the  other  three  effects  previously  mentioned,  and  partly 
because  the  theoretical  predictions  that  must  precede  an  experimental  test 
have  not  been  adequately  developed. 

Because  the  four  clearly  distinct  phenomena  can  cause  a  multiple  on¬ 
set  of  ground  backscatter,  investigators  should  approach  with  caution  the 
interpretation  of  a  measurement  of  the  leading  edge  of  ground  scatter. 

If  data  from  only  a  fixed  frequency  are  available,  it  seems  doubtful  (to 
this  author,  at  least)  that  one  could  discriminate  between  the  possible 
causes.  Such  discrimination  may  be  possible  with  sweep-frequency  data 
because  of  differing  systematic  variations  with  change  of  frequency.  For 
example,  the  magnetoionic  effect  and  the  AT-AR  effect  must  decrease  with 
an  increase  in  frequency,  but  the  synthesis  of  Fig.  8d  shows  that  the  lay- 
erlet  effect  increases  with  frequency. 

The  concentric  anomaly  is  unrealistic  because  it  covers  the  entire 
path.  It  seems  more  likely  that  such  an  anomaly  would  be  of  limited  hori¬ 
zontal  extent.  One  may  examine  the  rays  computed  from  the  concentric  m 'ri¬ 
els  to  deduce  the  effect  of  a  limited  anomaly,  noting  that  rays  that  do 
not  pass  through  the  anomaly  cannot  be  influenced  by  it.  Application  of 
this  simple  logic  shows  that  a  "concentric  anomaly"  of  limited  horizontal 
extent  would  create  an  extra  backscatter  leading  edge  over  a  limited  fre¬ 
quency  interval.  The  ray  calculations  show  that,  at  lower  frequencies, 
the  rays  that  are  able  to  penetrate  to  the  level  of  the  anomaly  leave  be¬ 
fore  reaching  the  range  of  the  anomaly.  When  the  frequency  is  too  high, 
the  rays  that  reach  the  anomaly  are  not  refracted  back  to  the  earth. 

The  experimental  data  show  that  the  false  leading  edges  usually  exist 
over  only  a  limited  interval  of  frequencies,  unlike  the  synthetic  examples. 
It  is  possible  that  most  of  the  experimental  measurements  of  multiple  lead¬ 
ing  edges  are  attributable  to  bounded  anomalies  which  are  a  little  more 
than  discontinuities  in  the  electron  density  vs  height  profile.  These  may 
or  may  not  be  moving  vertically.  Repetitive  backscatter  records  would  show 
vertical  motion  if  it  existed,  but  unfortunately,  the  experimental  data 
previously  mentioned  do  not  contain  any  multiple  leading  edges. 
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V  DRAWINGS  OF  THE  FOCUSED  RAYS 


Once  rays  are  calculated,  it  is  comparatively  easy  to  have  them 
automatically  drawn  by  the  plotting  equipment  that  is  associated  with 
most  large  computers.  These  plots  are  surprisingly  useful  because  they 
convey  much  information  in  a  compact,  readily  understood  way.  Three  ex¬ 
amples  showing  18-MHz  rays  are  presented  on  Fig.  9.  Part  (a)  was  com¬ 
puted  with  the  ambient  model,  and  parts  (b)  and  (c)  with  the  outward 
tilt  model.  In  all  three,  the  bottom  line  represents  the  surface  of  t!  a 
earth  with  range  marked  in  100-km  increments.  Rays  originate  at  the 
left  side,  leaving  the  earth  at  successive  half  degrees  of  takeoff  angle. 
On  each  ray  is  placed  a  tick  mark  orthogonal  to  the  ray  at  every  loca¬ 
tion  where  the  accumulated  group  delay  is  a  multiple  of  100  psec.  From 
inspection  of  the  pattern  of  the  ticks,  one  can  perceive  surfaces  of 
constant  delay,  which  might  be  considered  to  indicate  the  successive 
positions  of  a  short  pulse  at  100-psec  intervals.  (This  plot  can  be 
made  by  use  of  phase  instead  of  group  delay,  in  which  case  the  surfaces 
represent  wavefronts.) 

Signal  strength  can  be  deduced  from  these  drawings  if  one  visualizes 
the  rays  as  flux  lines  in  an  energy  flow.  This  concept  must  be  used  with 
moderate  caution  because  it  is  often  incorrect;  for  example,  neighboring 
rays  frequently  cross  and  ray  theory  does  not  work  in  such  a  circumstance. 
Nevertheless,  the  flux  line  idea  works  in  such  a  wide  region  that  it  is  a 
useful  concept. 

The  rays  in  the  ambient,  Fig.  9a,  show  that  the  energy  strikes  the 
distant  ground  uniformly  in  both  space  and  time  because  both  the  ray  spac- 
ings  and  the  delay  surface  spacings  are  uniform.  When  the  outward  tilt 
is  added,  spatial  focusing  occurs,  as  shown  in  Fig.  9b.  The  rays  converge 
as  they  strike  the  ground  near  2550  km,  causing  strong  backscatter  to  be 
detected  from  that  range.  Those  rays  that  take  off  at  low  angles  and  pass 
through  the  anomaly  on  their  upward  transit  strike  the  ionosphere  beyond 
the  anomaly  at  an  angle  that  is  steeper  than  it  would  have  been  in  the 
absence  of  an  anomaly  because  the  rays  are  refracted  downward  less  while 
they  are  in  the  anomaly  since  the  electron  density  there  is  lower  than  am¬ 
bient.  Since  these  are  "lower  rays,"  steeper  incidence  produces  shorter 
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delay  increments  of  100  j^sec  to  show  the  time  distribution  of 


range  and  it  can  be  seen  that  the  low  angle  rays  go  a  lesser  distance 
because  of  the  anomaly.  The  rays  that  encounter  the  anomaly  at  apogee 
are  not  much  affected  by  it.  They  still  go  about  2500  to  2600  km,  as 
they  did  in  the  ambient  case,  although  the  ray  spacing  is  disturbed. 
However,  the  short-range  rays  that  strike  the  anomaly  on  their  downward 
leg  are  not  refracted  downward  as  strongly;  consequently,  they  go  farther 
than  do  their  counterparts  in  the  ambient  ionosphere. 

The  combination  of  rays  appears  to  be  focused  because  the  long  dis¬ 
tance  rays  are  shortened  and  the  short  distance  rays  are  lengthened  in 
this  manner.  To  illustrate  the  mechanism  more  clearly,  Fig.  9  shows  the 
calculation  performed  with  a  90  percent  depletion  in  the  outward  tilt. 
This  exaggerates  the  effect  to  an  unrealistic  degree  but  serves  to  em¬ 
phasize  the  distortion  in  the  ray  pattern.  A  few  of  the  higher-angle 
rays  have  been  left  off  this  plot  to  clarify  the  display. 
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VI  ADDITIONAL  CALCULATIONS  RELATED  TO  THIS  STUDY 


A  shortened  version  of  the  preceding  chapters  has  been  submitted  to 
the  Journal  of  Atmospheric  and  Terrestrial  Physics,  which  will  devote  a 
special  issue  to  a  symposium  at  which  this  paper  was  presented.  (See 
Acknowledgment  for  details.)  The  remainder  of  this  report  will  be  de¬ 
voted  to  computational  details  and  additional  calculations  that  will  be 
primarily  of  interest  to  specialists  who  may  wish  to  carry  out  further 
work  of  a  similar  nature.  Those  readers  who  are  interested  only  in  the 
physical  significance  of  the  results  may  choose  to  skim  lightly  through 
the  remaining  pages. 

A.  Other  Ionospheric  Models 

Only  five  ionospheric  models  were  previously  mentioned,  but  many 
more  were  used  during  various  phases  of  this  study.  Each  of  these  models 
was  assigned  a  three-digit  ionosphere  identification  (IID)  number.  This 
identification  system  will  be  introduced  here  because  it  will  be  useful 
ns  an  aid  in  tying  together  the  remainder  of  the  data.  The  ambient  model 
was  described  in  Chapter  III;  it  has  IID  number  165.  The  local  anomaly 
was  IID  329  and  there  was  also  a  more  exaggerated  anomaly  of  identical 
shape  in  which  the  central  density  was  only  10  percent  of  ambient;  this 
was  IID  254.  The  outward  tilt  model  was  IID  327,  and  an  exaggerated  ver¬ 
sion  of  it  (10  percent  central  density)  was  IID  325.  Similarly,  the  in¬ 
ward  tilt  was  IID  328;  its  exaggerated  version  was  IID  326. 

There  were  four  concentric  models,  one  of  which  had  its  base  at  an 
altitude  of  250  km,  as  previously  described.  This  was  IID  320.  A  second 
version,  IID  319,  was  identical  except  that  its  base  altitude  was  200  km. 
There  were  two  variations  of  this  latter  model,  which  were  designed  to 
reveal  a  backscatter  detail  that  will  be  discussed  later.  The  two  vari¬ 
ations  were  IID  341  and  342.  To  fully  specify  a  concentric  model,  it  is 
only  necessary  to  present  a  plot  of  electron  density  vs  altitude.  Four 
such  plots  given  on  Fig  10  thus  fully  define  the  four  concentric  anom¬ 
alies.  Table  1  summarizes  all  these  IID  numbers. 

Chapter  III  provided  descriptions  of  IID  165,  329,  327,  328,  and  320, 
in  that  order.  The  rayplots  of  Fig.  9  were  calculated  from  models  165, 

29  SEL-6S-029 


«> /v’’ Yv  v S 'J  V,  v  V '  V  VJ /  ■.  a  f.  A.y\  ci Ki  r.  ci  vo  A  wni 1  • . 


■.  *-  *V  V 


300 

j - f - 1 - 1 - \ - - 

300|  "T“  T  ~\ — 

IID  341 

IID  342 

r  — 

1  1 

1  / 

1  /  I 

**0-  J 

|»  /  J 

E 

*10  — 

*io  — 

£  200—  _ 

200—  *00—  — 

V 

1 

h  4 

j-  /  X 

i 

P  -i 

r  /  i 

£ 

\  s'  1 

1  s  l 

P  ■] 

r  "1 

too*-  — 

100-  — * 

1 

i  ^ 

60 

- 1  <  ■  _ L_ _ 1 

60  1  _  1  I  1 

10*  10*  10*  10* 

io*  io4  W*  to* 

N.  Electron  Density  (per  cm*) 

N,  Electron  Density  (per  cm*) 

Fig.  10.  FOUR  CONCENTRIC  ANOMALY  MODELS. 


327,  and  325.  In  addition  to  the  models  outlined  in  Table  1,  the  local 
and  tilted  anomalies  were  tried  with  a  central  density  of  90  percent  of 
ambient.  These  models  did  not  produce  ray  disturbances  sufficiently 
strong  to  permit  the  easy  identification  of  anomalous  effects.  Appar¬ 
ently,  the  noise  in  the  computation  is  comparable  to  the  effect  of  the 
90  percent  anomalies  and  presumably  this  could  have  been  remedied  by  cal¬ 
culating  with  smaller  steps. 
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Table  1 


SUMMARY  OF  I ID  NUMBERS 


It  should  be  noticed  that  the  local  and  tilted  anomalies  are  all  re¬ 
gions  in  which  the  electron  density  is  more  rarified  than  in  the  ambient 
model.  In  general,  radio  rays  are  refracted  away  from  regions  of  elec¬ 
tron  concentration  and  thus  the  rays  would  be  drawn  toward  the  center  of 
a  rarified  region  and  it  would  act  as  a  converging  lens  acts  upon  light. 
The  idea  of  a  rarified  region  is  sketched  in  Fig.  11a.  In  part  (b)  of 
the  figure  is  shown  a  dense  region  in  which  one  would  expect  that  the  rays 
would  diverge.  This  was  not  simulated  in  this  study  and  so,  in  a  sense, 
we  ha\e  attacked  only  one  of  the  two  fundamental  types  of  anomaly.  The 
limitation  does  not  apply  to  the  concentric  models,  which  contain  mixtures 
of  over-  and  under-dense  regions. 
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RARIFIED  REGION 


a.  Rarified  anomaly  in  which  passing  rays  are  refracted  toward 
the  center 


b.  Dense  anomaly  in  which  rays  refract  away  from  the  center 

Fig.  11.  TWO  TYPES  OF  ANOMALIES. 

B.  Drawings  of  Rays  in  Various  Anomalies 

One  of  the  most  useful  forms  of  computer  output  is  the  ray  drawing, 
examples  of  which  have  been  presented  in  Fig.  9.  The  main  advantage  of 
this  form  of  output  apparently  stems  from  the  fact  that  the  data  format 
is  intuitively  obvious.  The  main  limitation  of  this  form  of  data  is  that 
it  shows  effects  for  only  a  single  frequency  on  a  single  plot.  The  ray 
drawings  also  suffer  from  the  disadvantage  that  they  cannot  show  very 
9 mall  perturbations  in  pattern  uniformity  because  the  drawing  is  neces¬ 
sarily  analog  and  the  scale  is  quite  large. 
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To  show  frequency  effects,  one  can  show  a  number  of  ray  drawings  in 
the  same  ionosphere  at  different  frequencies.  This  technique  is  used  on 
Fig.  12,  which  shows  rays  at  5  frequencies  in  IID  329,  the  local  anomaly 
with  a  central  density  70  percent  of  ambient.  In  the  third  and  fourth 
plots  on  this  drawing,  it  will  be  noticed  that  the  rays  do  not  emanate 
from  a  single  point  and  that  they  do  not  return  to  the  earth  but  rather 
miss  it  by  a  few  kilometers.  These  are  two  manifestations  of  an  error 
in  placement  of  the  plotting  pen  on  the  paper.  The  numerical  data  com¬ 
puted  for  these  rays  were  unaffected  by  this  source  of  error.  The  first 
and  third  plots  on  the  page  lack  any  upper  rays  but,  again,  this  is  not 
a  significant  deficiency  because  the  other  rays  were  calculated  at  a  dif¬ 
ferent  time  and  thus  they  could  not  be  put  on  the  same  plot. 

One  can  more  readily  see  the  significant  effects  on  those  plots  by 
placing  the  eye  near  the  plane  of  the  paper  on  the  right  side  and  looking 
back  lengthwise  along  the  plot.  This  has  the  effect  of  compressing  the 
horizontal  axis  and  emphasizes  vertical  discontinuities.  For  oblique 
rays  such  as  this,  most  of  the  refraction  effects  show  up  as  vortical 
changes.  In  addition,  however,  there  are  changes  in  the  time  distribu¬ 
tion  of  energy  which  can  be  seen  in  the  pattern  of  tick  marks.  For  exam¬ 
ple,  at  the  discontinuity  in  the  ray  pattern,  there  is  also  a  transverse 
shift  in  the  time  mark  pattern,  indicating  that  the  energy  that  impinges 
on  the  ground  is  not  only  focused  spatially  but  is  also  redistributed  in 
time.  Since  the  backscatter  measurements  are  influenced  by  both  temporal 
and  spatial  variations,  the  interruption  in  the  smooth  tick  mark  patterns 
is  significant  in  an  analysis  of  backscatter. 

The  top  drawing  on  Fig.  12  was  made  with  the  local  anomaly  located 
at  a  distance  of  only  1000  km.  The  remaining  four  ray  drawings  had  the 
anomaly  located  at  1200  km.  The  region  of  maximum  ray  pattern  discontin¬ 
uity  is  roughly  twice  the  distance  of  the  anomaly  itself,  except  in  the 
26-MHz  example  Where  the  skip  distance  exceeded  2  x  1200  km. 

Figure  13  shows  rays  in  the  exaggerated  local  anomaly,  IID  254.  In 
all  four  cases  illustrated,  the  anomaly  was  only  1000  km  from  the  trans¬ 
mitter.  By  again  sighting  lengthwise  down  the  drawings,  it  is  easier  to 
see  the  reason  for  the  specific  form  of  the  ray  discontinuities.  There 
appears  to  be  a  considerable  randomness  in  the  ray  pattern  in  this  set  of 
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Fig.  12.  RAYS  IN  THE  LOCAL  ANOMALY. 


calculations.  These  rays  were  traced  by  a  computer  program,  which  makes 
an  approximation  that  does  not  work  well  when  the  electron  density  is  high 
and  its  gradient  is  also  large  and  nearly  horizontal,  The  10  percent  anom¬ 
aly  models  (I  ID  254,  325,  and  326)  probably  represent  a  case  in  which  this 
raytracing  program  has  marginal  usefulness  because  these  models  possess 
severe  horizontal  gradients  that  may  account  for  the  disorder  in  these 
rays.  Nevertheless,  by  examining  them,  one  can  clearly  see  the  tendency 
of  the  rays  to  refract  toward  the  centerline  of  the  anomaly  which  has  its 
center  at  an  altitude  of  200  km  and  at  a  range  of  1000  km. 

To  understand  this  refraction  effect,  one  must  appreciate  the  differ¬ 
ence  between  upper  and  lower  rays.  These  two  types  of  rays  are  usually 
not  well  defined  but  the  author  has  found  the  following  definition  to  be 
useful:  If  a  ray  in  a  given  layer  would  have  a  shorter  range  if  its  take¬ 

off  angle  were  increased,  then  that  ray  is  a  lower  ray.  If  the  hypotheti¬ 
cal  increase  in  takeoff  angle  produced  an  increase  in  range,  then  the  ray 
is  an  upper  ray.  Stated  mathematically,  a  lower  ray  has  a  negative  dR/d|3 
while  an  upper  ray  has  a  positive  dR/d0. 

To  apply  this  upper  ray-lower  ray  reasoning  to  the  drawings  of  Fig. 

13,  notice  that  a  ray  that  is  traveling  along  an  upward  leg  of  an  iono¬ 
spheric  transit  may  be  refracted  either  upward  or  downward  by  the  local 
anomaly,  depending  on  whether  it  passes  above  or  below  the  anomaly  center¬ 
line.  If  the  ray  is  refracted  upward,  then  it  enters  a  trajectory  that 
corresponds  to  a  higher  takeoff  angle.  Therefore,  if  a  ray  is  a  lower  ray 
and  is  refracted  upward,  its  range  will  be  shortened.  Thus  we  see  that 
there  are  two  binary  criteria  and  consequently  there  are  four  possible 
cases : 

(1)  A  lower  ray  refracted  upward  will  have  decreased  range. 

(2)  A  lower  ray  refracted  downward  will  have  increased  range. 

(3)  An  upper  ray  refracted  downward  will  have  decreased  range. 

(4)  An  upper  ray  refracted  upward  will  have  increased  range. 

These  relations  usually  hold  true,  but  in  some  cases  the  refraction 
may  actually  change  a  ray  from  upper  to  lower.  As  a  consequence,  the  rea- 
soning  becomes  comparatively  complex  when  the  anomaly  under  study  sits  at 
half  the  skip  distance,  because  then  the  rays  that  pass  through  the  anom- 
aiy  are  readily  changed  from  upper  to  lower  (or  vice  versa)  by  a  slight 
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refraction.  Early  attempts  by  the  author  to  calculate  rays  in  local  anom¬ 
alies  led  to  some  confusion  which  may  now,  with  hindsight,  be  partly  at¬ 
tributed  to  positioning  of  the  anomaly  at  half  the  skip  distance.  This 
destroys  the  simple  patterns  of  ray  disturbance  that  have  been  shown  in 
the  preceding  calculations.  It  may  be  that  the  anomaly  effects  on  rays 
really  do  disappear  at  short  ranges,  for  examination  of  the  backscatter 
records  included  in  this  report  will  show  that  there  are  not  many  streaks 
at  short  ranges.  Most  are  beyond  600  km. 

With  the  above  thoughts  in  mind,  examine  the  20-MHz  rays  on  Fig.  13. 
Notice  that  the  locus  of  apogees  appears  to  follow  three  separate  arcs. 

The  first  arc  has  its  peak  at  a  range  of  about  900  km,  the  ^cond  arc 
peaks  at  about  1250  km,  and  the  third  arc  peaks  at  about  1600  km.  These 
three  arcs  intersect  at  roughly  1100  and  1500  km.  The  farthest  arc  is 
composed  of  those  rays  that  miss  the  anomaly  altogether,  having  passed 
underneath.  The  middle  arc  is  composed  of  the  apogees  of  those  rays  that 
passed  through  the  bottom  half  of  the  anomaly.  Notice  the  pattern  where 
these  rays  strike  the  ground.  There  is  a  deficiency  of  energy  from  about 
2300  to  about  2700  km  where  the  rays  diverge  from  one  another.  The  cause 
of  this  divergence  can  be  traced  backward  to  the  anomaly  where  it  is  seen 
that  the  rays  are  being  drawn  inward  toward  the  center  of  the  anomaly, 
that  is,  refracted  upward.  Thus,  as  the  original  takeoff  angle  at  the 
starting  r»oint  gradually  increases,  the  range  gradually  decreases  until 
the  rays  begin  to  strike  the  anomaly.  Then  the  effective  takeoff  angle 
increases  still  more  rapidly  so  that  the  ground  range  decreases  still 
more  rapidly  and  this  leads  to  the  deficiency  of  energy  seen  in  the  plot. 

The  shortest— range  arc  of  apogees  is  formed  of  those  rays  that  pass 
through  the  top  of  the  anomaly.  The  pattern  of  these  rays  on  the  ground 
is  not  so  well  defined  because  they  strike  near  the  skip  distance  where 
slight  deviations  change  the  rays  from  upper  to  lower  or  vice  versa.  The 
local  refraction  of  the  rays  while  they  are  in  the  anomaly  can  clearly  be 
seen  in  the  upper  rays  at  18  and  20  MHz.  Because  the  local  anomaly  is  of 
parabolic  cross  section  (in  electron  density  vs  distance),  the  highest 
electron  density  gradient  occurs  at  the  outer  edge  of  the  anomaly;  this 
is  where  the  inward  refraction  will  be  highest.  Thus  the  downward  refrac¬ 
tion  of  the  upper  rays  is  clearly  seen  to  occur  at  an  altitude  just  short 
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Fig.  8b.  The  18-  MHz  example  was  also  shown  as  part  (b)  of  Fig.  9.  From 
this  family  of  rays,  one  can  see  the  tendency  for  the  range  of  the  focus¬ 
ing  legion  to  increase  with  increasi)  g  frequency.  It  is  also  possible  to 
gain  some  insight  about  the  focusing  mechanism,  although  this  is  more  eas¬ 
ily  done  with  the  ray  drawings  of  Fig.  15  computed  in  the  exaggerated 
equivalent  of  the  same  ionospheric  model.  Here,  fewer  rays  are  shown 
and  also  the  upper  rays  are  not  included.  The  18-MHz  example  is  the  same 
as  that  which  was  shown  as  part  (c)  of  Fig.  9.  By  sighting  lengthwise 
along  the  line  of  apogees,  one  can  see  an  upward  bulge  in  each  ray  family 
just  beyond  the  region  where  the  outward  tilted  anomaly  passes  through  the 
apogees.  Using  the  upper  ray-lower  ray  guidelines  given  above,  one  can  de¬ 
duce  the  mechanism  of  focusing.  This  was  described  at  length  during  the 
discussion  of  Fig.  9,  but  it  is  perhaps  more  clearly  illustrated  here 
where  several  examples  can  be  compared.  The  basic  reason  for  the  upward 
bulge  in  the  line  of  apogees  is  the  fact  that  the  electron  density  de¬ 
creased  so  that  in  this  region  there  was  less  barrier  to  the  penetration 
of  rays. 

Figure  16  shows  a  number  of  ray  calculations  in  the  inward  tilt  model, 
I  ID  328.  This  should  be  compared  to  Fig.  14  to  see  the  decrease  in  range 
of  ray-focusing  vs  frequency  as  contrasted  to  the  increase  when  the  anom¬ 
aly  has  an  outward  tilt.  The  exaggerated  anomaly,  IID  326,  was  used  to 
generate  the  rays  given  on  Fig.  17,  which  should  be  compared  with  Fig  15. 

By  again  examining  the  line  of  apogees,  one  can  see  the  qualitative  dif¬ 
ference  between  the  rays  in  the  inward  tilt  and  those  in  the  outward  tilt. 
The  inward  tilted  anomaly  is  inclined  toward  the  rays  and  so  the  beginning 
of  its  effect  tends  to  be  more  abrupt  while  its  ending  tends  to  be  more 
smooth  than  the  corresponding  features  in  the  apogee  lines  on  Fag.  15. 

To  see  this,  notice  that  the  upward  bulge  in  the  apogee  lines  on  Fig.  15 
is  almost  symmetric  with  respect  to  the  bulge  centerline.  In  Fig.  17,  the 
upward  bulge  has  a  comparatively  sharp  left  side,  but  the  right  side  is  so 
gradual  that  it  is  not  clearly  defined.  Except  for  this  minor  difference, 
the  focusing  mechanism  for  the  two  kinds  of  tilted  anomalies  appears  to  be 
essentially  identical. 

Figure  18  shows  an  interesting  comparison  of  rays  at  the  same  frequen¬ 
cy  in  several  different  models.  As  a  reference  standard,  the  top  figure 
shows  ravs  in  the  ambient  ionosphere  and  the  upper  rays  have  been  deleted 
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Fig.  18.  RAYS  AT  18  MHz  IN  MANY  DIFFERENT  IONOSPHERIC  MODELS. 


in  order  to  emphasize  the  uniformity  in  the  lower  ray  family.  It  is  pos¬ 
sible  to  detect  some  slight  jitter  in  these  rays;  all  such  jitter  is  due 
to  computation  defects.  The  programmed  calculation  was  identical  to  that 
used  in  producing  the  other  ray  drawings  of  this  report. 

Comparison  of  the  lines  of  apogees  for  the  three  models  with  a  10 
percent  central  density  will  serve  to  emphasize  the  three  different  quali¬ 
tative  behaviors  of  the  rays  in  the  three  different  anomaly  types.  For 
the  local  anomaly,  the  line  of  apogees  is  formed  of  three  arcs.  With  the 
outward  tilt,  the  line  of  apogees  has  a  single  upward  hump,  which  is  rough¬ 
ly  symmetric  about  its  own  centerline.  For  the  inward  tilt,  the  hump  in 

the  line  of  apogees  is  sharply  defined  on  the  left  side  but  less  well  de¬ 

fined  on  the  right. 

C.  Backscatter  under  Various  Circumstances 

The  set  of  backscatter  calculations  presented  in  Chapter  IV  had  a 
logarithmic  frequency  scale  and  a  dynamic  range  of  30  dB.  Many  other 
calculations  have  been  carried  out;  some  of  the  more  useful  examples 
will  be  presented  here.  There  are  three  additional  classes  of  back¬ 
scatter  records: 

(1)  backscatter  from  other  ionospheric  models  besides  those  used 
in  Chapter  IV, 

(2)  calculations  with  the  dynamic  range  set  at  60  dB,  and 

(3)  display  with  a  linear  frequency  scale  that  changes  the  apparent 

echo  shape. 

Figure  19  shows  backscatter  from  the  exaggerated  versions  of  the 
first  three  anomalies  used  to  generate  Fig.  8.  All  the  calculation  pa¬ 
rameters  are  the  same  and  the  figure  layout  is  similar,  to  facilitate 
comparison,  but  calculations  in  the  exaggerated  anomalies  were  carried 
out  only  at  frequencies  of  12,  14,  16,  18,  and  20  MHz.  (The  program  that 
plots  the  backscatter  automatically  adjusts  the  amplitude  of  each  single¬ 
frequency  record  so  that  it  has  an  abscissa  span  equal  to  that  of  a  sin¬ 
gle  megahertz.)  It  can  be  seen  that  the  exaggerated  anomaly  creates  a 
very  strong  disturbance  in  the  backscatter,  so  strong  in  fact  that  there 
is  no  longer  a  clearly  identifiable  streak  present. 

Figure  20  shows  the  backscatter  obtained  with  three  different  con¬ 
centric  anomalies.  Part  (a)  of  the  figure  shows  the  result  when  the 
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Fig.  19.  BACKSCATTER  IN  THE  EXAGGERATED  ANOMALIES, 
(a)  With  the  exaggera  ed  local  anomaly;  (b)  with 
the  exaggerated  outward  tilt;  and  (c)  with  the 
exaggerated  inward  tilt. 


full-cycle  sinusoid  density  multiplier  of  Fig.  7a  is  moved  downward  50  km 
relative  to  the  location  which  was  previously  described.  The  backscatter 
computed  with  this  model  is  shown  in  Fig.  20a  and  should  be  compared  to 
that  shown  in  Fig.  8d,  which  is  otherwise  the  same.  Notice  that  there  are 
two  major  new  effects  on  Fig.  20a,  not  on  Fig.  8d: 

(l)  The  false  leading  edge  appears  at  a  lower  frequency,  but  the 

reduction  in  minimum  time  delay  is  not  nearly  as  marked  at  the 
high  frequencies.  In  fact,  the  effect  of  the  anomaly  practi¬ 
cally  disappears  above  25  MHz  in  Fig  20a,  whereas  the  effect 
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Fig.  20.  BACKSCATTER  IN  OTHER  CONCENTRIC  ANOMALIES. 

(a)  Full  sinusoid  multiplier  with  base  at  200  km; 

(b)  positive  half  cycle  of  sinusoid,  electrons  added 
above  210  km;  and  (c)  negative  half  cycle  of  sinus¬ 
oid,  electrons  deleted  above  200  km. 

of  the  anomaly  exerted  dominant  control  over  the  backscatter  at 
high  frequencies  in  Fig.  8d .  This  different  relative  behavior 
is  attributable  to  the  fact  that  the  higher  frequency  rays  have 
higher  apogees,  which  are  more  strongly  influenced  by  the  higher 
anomaly . 

(2)  There  is  another  apparent  leading  edge  within  the  body  of  the 
backscatter,  which  first  appears  at  approximately  13  MHz  just 
beyond  the  strongest  peak  in  the  backscatter.  At  higher  fre¬ 
quencies,  the  effect  becomes  more  pronounced,  and  its  time  delay 
is  more  rapidly  increasing  with  respect  to  frequency  than  that 
of  either  of  the  other  two  signal  enhancements.  The  effect 
finally  disappears  at  23  MHz. 
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In  an  effort  to  understand  why  these  two  separate  effects  should  oc¬ 
cur,  the  raysets  were  examined  and  the  pattern  of  the  numbers  appeared 
to  indicate  that  the  first  leading  edge  could  be  attributed  solely  to 
the  positive  half  cycle  of  the  sinusoid,  while  the  last  leading  edge 
could  be  attributed  to  the  negative  half  cycle  (that  is,  to  the  rare¬ 
faction).  It  is  difficult  to  be  certain  of  this  cause-effect  relation 
from  an  examination  of  tabulated  numbers,  and  so  a  chock  was  conducted 
by  making  two  ionospheric  models,  one  with  the  positive  hnlf  cycle  and 
the  other  with  a  negative  half  cycle.  These  are  I  ID  341  and  342,  which 
have  previously  been  described. 

With  the  positive  half  cycle  (electrons  added  between  210  and  220  km) , 
we  obtained  the  backseatter  shown  on  Fig.  20b.  With  the  negative  half 
cycle  (electrons  removed  between  200  and  210  km) ,  we  obtained  the  back- 
scatter  of  Fig.  20c.  Comparison  of  these  two  backseatter  plots  with 
that  of  Fig.  20a  shows  that  there  is  indeed  a  tendency  for  the  positive 
half  cycle  to  produce  the  added  leading  edge  at  shortest  time  delay  and 
for  the  negative  half  cycle  to  produce  only  the  longest  time  delay  effect. 
However,  the  separation  of  the  two  effects  is  not  complete.  In  Fig.  20b, 
there  is  some  evidence  of  a  third  leading  edge,  which  begins  at  19  MHz 
and  extends  up  through  24  MHz,  a  higher  frequency  than  before.  In  Fig. 

20c,  the  fully  developed  effect  is  seen  just  as  it  was  in  part  (a)  of 
the  figure,  but  in  addition  there  is  a  new  signal  enhancement,  which  ap¬ 
pears  at  19  MHz  and  can  clearly  be  identified  through  24  MHz. 

In  some  of  the  charts  that  follow,  raypath  parameter  relations  will 
be  presented  which  will  clarify  the  mechanisms  underlying  these  baek- 
seatter  discontinuities.  At  this  point,  we  have  only  seen  that  if  wo 
insert  a  particular  ionospheric  model  into  the  computer,  out  comes  a  new 
and  different  backseatter  anomaly. 

Backseatter  with  60-dB  Dynamic  Range.  In  Chapter  II,  Part  C,  it  was 
explained  that  the  backseatter  simulation  program  has  a  variable  "dynamic 
range,"  which,  in  fact,  controls  the  total  number  of  decibels  of  the  com¬ 
puted  backseatter  that  will  be  presented  on  the  display.  After  the 
amplitude  is  computed  at  each  fixed  frequency,  all  other  amplitudes  arc 
compared  with  the  maximum.  If  the  dynamic  range  setting  is  30  dB,  for 
example,  then  r.he  computer  is  instructed  to  disenrd  any  echo  that  is 
not  within  30  dB  of  the  maximum.  For  example,  if  the  peak  backseatter 
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echo  power  were  10  W,  then  the  display  would  present  only  echo  compo- 

-10  -7 
nents  that  exceed  10  W,  that  is,  30  dB  less  than  10 

All  of  the  preceding  data  has  been  calculated  with  a  dynamic  range 
setting  of  30  dB.  In  this  section,  we  will  show  many  of  the  same  records 
calculated  with  a  60-dB  setting.  This  tends  to  produce  unrealistic  ef¬ 
fects  because  the  approximate  backscatter  simulation  technique  introduces 
larger  errors  than  does  the  more  rigorous  technique  described  in  Ref.  1. 

Nevertheless,  the  60-dB  setting  can  be  useful  because  some  of  the  echo 
components  lie  between  30  and  60  dB.  The  sharp  trailing  edges  seen  in 
the  60-dB  backscatter  are  not  real  but  are  caused  by  approximations.  In 
most  of  this  backscatter,  the  amplitude  should  die  off  gradually  with  in¬ 
creasing  time  delay,  much  as  it  appears  to  do  on  the  30-dB  records.  By 
computing  the  original  rays  with  a  finer  integration  step  and  using  finer 
time  steps  in  the  backscatter  synthesis,  it  should  be  possible  to  achieve 
better  realism  in  the  60-dB  records  by  use  of  these  same  programs  without 
any  modification.  However,  we  have  not  attempted  to  do  this  since  it  is 
our  general  practice  to  use  the  more  rigorous  backscatter  simulation 
technique  for  the  examination  of  smaller  details. 

Figure  21a  shows  60-dB  backscatter  in  the  ambient  model  and  can  best 
be  compared  to  the  example  shown  in  Fig.  5.  At  the  higher  frequencies  on 
the  60-dB  model,  it  is  seen  that  the  trailing  edge  has  the  artificial 
sharp  cutoff  previously  mentioned.  Also,  at  the  low  signal  levels  at  the 
lower  frequencies,  it  is  clear  that  the  stairstep  structure  in  the  echo 
amplitude  vs  time  delay  is  attributable  to  the  finite  number  of  rays  that 
were  calculated.  Neither  of  these  effects  would  appear  if  we  had  calcu¬ 
lated  a  very  large  number  of  rays  originally.  It  can  be  seen  from  Fig.  5 
that  the  calculation  parameters  were  adequate  for  a  30-dB  presentation, 
however. 

Parts  (b) ,  (c)  and  (d)  of  Fig.  21  should  be  compared  with  Fig.  8. 

With  a  higher  dynamic  range,  it  is  possible  to  see  the  streak  structures 
over  a  wider  frequency  range.  For  example,  in  the  outward  and  inward 

tilts,  the  streak  is  present  at  8  MHz  but  it  is  so  weak  that  it  can  be  * 

seen  only  on  the  records  with  a  higher  dynamic  range.  Similarly,  the 
local  anomaly  creates  a  weak  effect  at  12  MHz,  which  is  seen  only  on  the 
60-dB  record. 
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Fig.  21.  BACKSCATTER  IN  AMBIENT  AND  LOCALIZED 
ANOMALIES  WITH  60  dB  DISPLAYED.  (a)  With  the 
ambient  ionospheric  model;  (b)  with  the  local 
anomaly;  (c)  with  the  outward  tilt;  and 
(d)  with  the  inward  tilt. 
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Fig.  22.  BACKSCATTER  WITH  CONCENTRIC  IONOSPHERIC 
MODELS  DISPLAYED  TO  60  dB.  (a)  With  the  full  cycle 
sinusoid  based  above  250  km  and  (b)  with  the  full 
cycle  sinusoid  based  above  200  km. 

Figure  22  shows  60-dB  backscatter  in  the  two  concentric  ionospheres, 
one  with  the  anomaly  based  at  250  km  and  the  other  with  the  anom. based 
at  200  km.  Little  additional  information  is  seen  on  these  records;  they 
serve  but  to  show  that  the  30-dB  dynamic  range  setting  was  a  better  choice. 
In  fact,  on  part  (b)  of  Fig.  22,  the  third  signal  onset  is  just  barely  vis¬ 
ible  whereas  the  same  third  onset  is  clearly  seen  in  the  30-dB  record  shown 
on  Fig.  20a.  Thus,  in  this  case,  the  30-dB  setting  shows  all  the  useful 
information  more  clearly  than  does  the  60-dB  recon1.  This  is  an  interest¬ 
ing  example  which  shows  that  too  much  dynamic  range  can  bo  detrimental. 

Figure  23  introduces  the  linear  scale  by  showing  the  srune  data  with 
both  a  logarithmic  and  a  linear  frequency  axis.  As  a  general  rule,  those 
who  do  analysis  usually  desire  to  have  a  logarithmic  frequency  axis  on  HF 
sounder  data  because  it  greatly  simplifies  many  calculations.  The  axis 
can  be  used  in  a  more  gel  aral  way  when  it  is  logarithmic  because,  in  a 
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Fig.  23.  COMPARISON  OF  BACKSCATTER  OBTAINED  WITH 
THE  AMBIENT  MODEL  DISPLAYED  WITH  LOGARITHMIC  AND 
LINEAR  FREQUENCY  SCALES.  (a)  With  a  logarithmic 
scale  and  (b)  with  a  linear  scale. 

sense,  it  operates  as  one  of  the  axes  of  a  slide  rule.  For  example,  if 
we  calculate  the  shape  of  sweep  frequency  backscatter  for  a  certain  iono¬ 
spheric  model  with  a  critical  frequency  of  10  MHz,  then  that  backscatter 
shape  is  identically  the  same  as  would  be  calculated  with  the  same  iono¬ 
sphere  if  all  electron  densities  were  cut  to  half  provided  that  all  radio 
frequencies  on  the  backscatter  record  are  cut  by  l/^/jf.  The  frequency 
scaling  can  be  done  with  a  simple  horizontal  shift,  provided  that  the 
frequency  scale  is  logarithmic.  This  scaling  works  because  we  neglect 
the  geomagnetic  field  so  that  the  index  of  refraction  is  simply 
Vl  -  80 . 6N/P ^ ,  and  thus  the  2’s  cancel. 

Those  whose  primary  occupation  is  the  gathering  of  sounder  data  tend 
to  iavor  the  use  of  a  linear  scale  because  i+  is  much  easier  to  construct 
equipment  that  will  produce  records  with  a  linear  scale.  As  a  consequence, 
most  experimental  backscatter  data  have  a  linear  frequency  scale.  This  is 

51  SFL-b8-0?9 


to) 

30f—  110  256 

20- 


lOg 


°L-^ 


l  I- 


10 


Rodio  Frequency  MHz 


unfortunate  for  analysis  because  the  shape  of  backscatter  depends  markedly 
on  whether  the  frequency  axis  is  linear  or  logarithmic.  If  the  electron 
density  everywhere  is  cut  in  half  and  wo  calculate  backscatter  with  a  lin¬ 
ear  frequency  scale,  then  not  only  do  the  frequencies  change  but  also  the 
overall  shape  of  the  backscatter  echo  will  change.  It  is  not  possible  to 
use  the  sliding  scale  concept  to  eliminate  the  added  complexity. 

In  order  to  serve  these  two  communities  of  users,  the  backscatter 
simulation  program  plots  synthetic  backsc.  i er  twice,  once  with  the  linear 
scale  and  once  with  the  logarithmic  sea  V  few  of  the  linear  examples 

are  shown  here  to  illustrate  the  changed  appearance.  On  Fig.  23b,  with 
the  linear  axis,  it  is  seen  that  t^e  peak  abscissa  displacement  of  each 
fixed  frequency  record  is  the  same  regardless  of  frequency,  because  now 
every  single  megahertz  spans  the  same  distance  on  the  abscissa.  In  sub¬ 
jective  appearance  this  is  an  advantage,  but  there  is  no  actual  practical 
benefit  due  to  this  new  equality  of  amplitudes  because  the  approximate 
backscatter  simulation  program  does  not  calculate  absolute  amplitude  in 
the  first  place.  A  normalization  has  already  taken  place  before  the  pro¬ 
gram  chooses  whether  to  use  a  linear  or  logarithmic  scale. 

Figure  24  shows  the  original  backscatter  set  from  Fig.  8  replotted 
with  linear  frequency  axes.  The  layout  is  identical  to  that  of  Fig.  8, 
but  no  new  information  is  added.  Notice  the  changed  appearance  of  the 
backscatter  with  the  linear  axis. 

D.  Reflectrix  Drawings  for  Different  Ionospheric  Models 

As  was  stated  previously,  the  ray  drawings  and  backscatter  drawings 
are  useful,  but  it  is  sometimes  difficult  to  understand  propagation  phe¬ 
nomena  simply  from  examining  these  particular  data  displays.  Another  kind 
of  data  display  that  often  lends  additional  insight  is  called  the  reflec¬ 
trix.  This  construction  was  introduced  by  Lejay  and  Lepechinsky  (Ref.  11) , 
who  coined  the  name  for  it.  In  fact,  the  reflectrix  has  been  defined  for 
use  only  with  a  symmetric  ionosphere  model  such  as  has  been  used  here  for 
the  ambient  ionosphere  and  for  the  concentric  anomalies.  With  the  local 
or  tilted  anomalies,  the  reflectrix  must  be  redefined. 

Figure  25  serves  to  define  the  reflectrix,  provided  the  electron 
density  model  is  concentric  with  the  earth.  When  the  ionosphere  contains 
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Fig.  24  BACKSCATTER  WITH  THE  ORIGINAL  SET  OF  ANOMALIES 
DISPLAYED  WITH  A  LINEAR  FREQUENCY  SCALE  (a)  With  a 
local  anomaly;  (b)  with  the  outward  tilt;  (c)  with  the 
inward  tilt  and  (d)  with  the  concentric  anomaly. 
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nonvertical  gradients  (tilts),  there  is  no  conventional  definition  and 
the  author  chooses  to  use  the  following:  "Range"  is  the  total  range  of 
the  ray,  and  "equivalent  height"  is  the  height  of  the  intersection  of  the 
two  straight  lines  which  are  extensions  of  the  ray  at  the  two  points  where 
it  encounters  the  earth.  In  this  case,  the  takeoff  angle  axis  is  neces¬ 
sarily  meaningless.  For  the  concentric  ionosphere,  it  was  possible  to 
have  all  three  axes  because  they  are  redundant;  that  is,  in  a  symmetric 
ionosphere,  if  one  knows  two  of  the  three  parameters  R,  H' ,  ft,  then  one 
can  calculate  the  third  without  further  information.  In  a  tilted  iono¬ 
sphere,  this  is  no  longer  true;  as  a  consequence,  one  can  no  longer  use 
all  three  axes  simultaneously.  This  is  why  the  definition  of  the  reflec- 
trix  has  been  restricted  to  concentric  ionospheres.  In  tilted  ionospheres 
one  could  choose  to  preserve  ft  and  R,  or  H'  and  ft,  or  H'  and  R.  We 
have  chosen  the  latter  alternative  here,  but  the  reader  is  cautioned  that 
this  is  not  a  universally  accepted  definition. 

Figure  26  shows  (a)  a  reflectrix  as  it  comes  from  the  automatic  plot¬ 
ter  and  (b)  a  version  that  has  been  cleaned  up  by  a  draftsman.  During  the 
automatic  plotting,  we  place  a  small  x  at  the  location  of  each  rayset 
and  then  draw  a  line  from  the  center  of  that  x  to  the  center  of  the  next 
x  which  is  located  by  the  next  rayset  in  a  deck  of  cards.  The  line  is 
drawn  only  when  the  two  successive  raysets  have  the  same  radio  frequency. 
The  line  that  joins  two  raysets  is  straight;  therefore  it  is  seen  in  the 
upper-ray  regions  where  the  rays  are  sparse  that  the  curve  is  no  longer 
apparently  smooth  but  is  composed  of  visibly  separate  straight-line  seg¬ 
ments.  Of  course,  this  is  an  artifice  that  can  be  eliminated  by  computing 
more  rays  in  the  upper  ray  region. 

The  automatically  drawn  plots  are  often  useful  because  one  can  ob¬ 
tain  a  measure  of  signal  strength  by  observing  the  closeness  of  spacing 
of  the  individual  x's.  In  the  drafting  process,  we  eliminate  the  x's 
and  substitute  a  smooth  curve,  which  is  easier  to  read,  but  it  lacks  an 
indication  of  the  takeoff  angle  spacing  of  the  rays  from  which  the 
reflectrix  was  constructed.  Because  this  section  of  the  report  is  in¬ 
tended  primarily  for  raytracing  specialists,  the  rcflectrixes  will  be 
presented  here  directly  as  they  were  plotted  without  any  redrafting, 
except  for  the  addition  of  a  few  labels  where  they  seem  advisable. 
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a.  The  reflectrixes  as  automatically  plotted 


b.  The  same  reflectrixes  drawn  by  a  draftsman 

Fig.  26.  REFLECTRIX  SHAPES  IN  THE  AMBIENT  IONOSPHERE. 

Both  reflectrixes  from  Fig.  26  were  calculated  for  the  ambient  model 
at  each  megahertz.  In  the  tilted  ionospheres,  it  was  our  practice  to  com¬ 
pute  only  the  even  megahertz  rays  as  an  economy  measure. 

Figure  27  shows  the  reflectrixes  in  the  local,  outward  tilt  and  in¬ 
ward  tilted  anomalies.  The  rays  for  frequencies  of  10  MHz  and  below  were 
not  calculated  in  the  local  anomaly  because  it  could  be  seen  in  the  am¬ 
bient  ionosphere  that  such  rays  would  not  have  reached  the  altitude  of 
the  anomaly  and  thus  would  not  have  been  affected  by  it.  From  the  com¬ 
parison  of  parts  (a) ,  (b) ,  and  (c)  of  Fig.  27  one  can  clearly  see  how  the 
range  and  time  delay  of  the  anomaly  disturbance  exhibit  a  frequency  de¬ 
pendence  that  is  markedly  influenced  by  the  shape  of  the  anomaly.  With  a 
local  anomaly,  part  (a),  the  disturbances  in  the  ray  pattern  are  located 
almost  directly  above  2400  km.  With  the  outward  tilt  in  part  (b)  of  the 
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a.  With  'he  local  anomaly 


b.  With  the  outward  tilt 


c.  With  the  inward  tilt 


Fig.  27  REFLECTRIXES  IN  THE  LOCALIZED  ANOMALIES. 


figure,  it  is  readily  apparent  that  the  effect  moves  to  a  greater  range 
at  a  greater  time  delay  with  an  increase  in  frequency  while  the  opposite 
trend  is  shown  with  the  inward  tilt  on  part  (c). 
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Figure  28  shows  the  reflectrix  calculations  for  the  exaggerated  ver¬ 
sions  of  the  local  and  tilted  anomalies.  It  is  probable  that  much  of  the 
fine  structure  on  these  curves  is  due  to  some  form  of  computation  error. 

In  particular,  the  erratic  behavior  of  the  rays  of  the  skip  distance  at 
a  few  selected  frequencies  is  probably  a  manifestation  of  some  form  of 
error  that  occurs  when  the  rays  pass  through  the  anomaly  near  their  apogee. 
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Notice,  for  example,  the  20-MIIz  ravs  in  Fig.  28a  and  28c.  It  seems  (on 
intuitive  grounds,  at  least)  that  the  disturbance  of  the  ray  pattern  should 
be  less  random  than  it  appears  in  these  examples.  The  computation  of  p 
ray  as  it  passes  through  an  anomaly  is  somewhat  difficult  because  the  pro¬ 
cess  is  necessarily  one  of  sampling  the  electron  density  and  its  gradient 
at  a  finite  number  of  locations  along  the  ravpath.  The  ray  structure  is 
influenced  by  the  way  in  which  t!  e  samples  are  s  ilected  and  it  appears 
that  we  have  not  found  an  optimum  solution  here.  The  errors  that  arise 
in  the  program  are  accentuated  in  the  exaggerated  disturbances,  which 
were,  in  fact,  constructed  primarily  to  test  the  severity  of  the  electron 
disturbances  that  could  be  accommodated  by  the  raytracing  program. 

Figure  29  shows  reflectrix  structures  in  the  four  concentric  models. 

Since  these  lack  any  horizontal  gradient,  the  original  reflectrix  defini¬ 

tion  of  Fig  25  is  fu'ly  applicable  and  the  takeoff  angle  axis  has  validity 
on  this  figure. 

The  effect  of  a  concentric  anomaly  is  a  small-scale  version  of  the 

effect  of  a  regular  ionospheric  layer.  To  see  what  happens  when  an  E 

layer  is  added  underneath  an  F  layer,  the  reader  can  refer  to  Ref.  12 
where  the  author  has  shown  a  variety  of  reflectrix  shapes.  Ideally,  the 
curve  which  characterizes  one  frequency  should  not  cross  any  other  curve 
which  applies  to  some  other  frequency.  We  see  in  Ref.  12  that  the  pres¬ 
ence  of  a  two-layer  structure  can  cause  just  such  a  crossing  to  occur. 

Thus,  this  undesirable  effect  occurs  in  a  circumstance  that  is  very  common. 
Many  other  reflectrix  plots  (not  shown  here)  were  almost  unreadable  be¬ 
cause  of  the  large  number  of  crossings  that  made  it  difficult  to  follow 
the  pattern  of  the  lines. 

E.  A  Modified  Form  of  the  Reflectrix 

Reference  to  the  reflectrix  explanation  of  Fig.  25  (or  to  the  more 
detailed  description  given  in  Ref.  12)  will  serve  to  show  that  the  dis¬ 
tance  from  the  origin  to  a  point  on  the  reflectrix  should  provide  a  mea¬ 
sure  proportional  to  the  group  delay  of  the  ray.  This  is  because  the 
theorem  of  Breit  and  Tuve  (Ref.  13)  states  that  the  group  delay  along  a 
ray  is  equal  to  the  delay  at  the  speed  of  light  along  the  equivalent 
straight-line  ray,  at  least  in  a  flat  earth-llat  ionosphere  situation. 
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a.  With  tht  full  cycle  sinusoid  based  above  250  km 
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b.  With  the  full  cycle  sinusoid  based  above  200  km 
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c.  Positive  half  cycle  of  sinusoid,  electrons  added  above  210  km 
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d  Negative  half  cycle  of  sinusoid,  electrons  deleted  above  200  km 
Fig  29  REFLECTRIXES  IN  THE  CONCENTRIC  ANOMALIES 
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This  is  approximately  true  in  a  curved  earth-curved  ionosphere  geometry, 
but  there  is  a  degree  of  error  that  is  worse  at  longer  distances. 

When  interest  is  centered  upon  ground  backscattcr,  then  one  of  the 
most  interesting  ray  parameters  is  the  range  vs  group  delay.  The  takeoff 
angle  has  a  comparatively  minor  significance,  and  an  approximate  measure 
of  takeoff  angle  would  be  sufficient.  In  the  reflectrix,  the  takeoff 
angle  information  is  exact,  but  the  time  delay  information  is  approximate. 
We  would  like  to  have  exact  time  delay  information  and  approximate  takeoff 
angle  information.  This  can  be  accomplished  by  a  slight  modification  of 
the  reflectrix  for  which  the  author  has  coined  the  term,  "timetrix."  In 
the  timetrix  construction,  the  distance  of  a  point  from  the  origin  of  the 
plot  is  fixed  by  the  distance  that  light  would  travel  in  half  the  transit 
time  of  the  ray.  For  this,  the  drawing  scale  of  Fig.  25,  part  (a),  (b), 
or  (c)  is  used.  Thus,  when  tnc  drawing  scale  is  doubled  as  it  is  in  the 
transition  to  part  (d)  of  Fig.  25,  then  the  distance  of  a  reflectrix  point 
from  the  origin  is  equal  to  the  total  distance  that  light  would  travel 
during  the  ray  transit  time,  that  is,  to  the  group  path  of  the  ray.  As 
in  the  reflectrix,  the  range  of  a  ray  is  used  as  the  second  parameter  fix¬ 
ing  the  location  of  a  point.  The  equivalent  height,  H * ,  is  now  slightly 
different  from  that  which  was  defined  in  Fig.  25,  the  difference  being  a 
direct  measure  of  the  error  in  the  theorem  of  Breit  and  Tuvc  when  it  is 
applied  to  curved  earth-curved  ionosphere  geometry.  Similarly,  there  is 
a  slight  error  in  the  takeoff  angle  measure  using  the  simple  protractor 
construction  of  Fig.  25d.  This  error  seldom  exceeds  3°,  however,  and  it 
is  this  severe  only  at  the  longer  ranges. 

Figure  30  shows  the  timetrix  for  the  ambient  ionosphere  and  should  be 
compared  with  Fig.  26,  which  shows  the  identical  data  in  reflectrix  con¬ 
struction.  On  Fig.  30,  at  any  frequency,  the  initial  backscattcr  will  be 
received  via  that  ray  that  is  nearest  the  origin;  this  can  be  located  by 
striking  an  arc  centered  on  the  origin  that  is  tangent  to  the  appropriate 
curve.  Analysis  of  backscattcr  can  be  aided  by  considering  the  intersec¬ 
tion  of  an  expanding  circle,  centered  on  the  origin,  with  one  of  these 
timetrix  curves.  One  can  thus  follow  the  progress  of  the  energy  by  all 
modes  as  n  function  of  time.  In  Ref.  12,  this  process  of  reasoning  was 
fully  explained  but  only  in  reference  to  a  reflectrix,  where  the  process 
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is  approximate.  The  interested  reader  may  refer  to  that  explanation 
which,  when  applied  to  the  timetrix,  is  exact. 

The  comparison  of  Fig.  30  with  Fig.  26b  can  be  used  to  establish  the 
amount  of  error  in  the  Breit  and  Tuve  theorem  in  the  following  manner: 
Select  any  point  and  locate  it  on  both  the  timetrix  and  the  reflectrix. 

■or  this  process,  it  is  easiest  to  use  a  high-frequency  upper-ray,  for 
there  the  error  is  greatest.  On  the  reflectrix,  measure  the  distance  of 
the  point  from  the  origin.  The  theorem  of  Breit  and  Tuve  would  tell  us 
that  this  is  the  time  delay  of  the  ray.  The  same  measurement  on  the  time¬ 
trix  yields  the  actual  time  delay  and  the  difference  of  the  two  measure¬ 
ments  is  then  the  error  in  the  theorem.  It  can  thus  be  seen  that  the 
theorem  is  quite  good  for  most  purposes,  particularly  when  one  considers 
how  little  is  known  about  the  ionospheric  structure  which  is  used  in  the 
original  generation  of  the  reflectrix  or  timetrix. 

Figure  31  shows  timetrix  structures  for  the  local  and  tilted  anoma¬ 
lies,  and  Fig.  32  shows  these  structures  for  the  concentric  anomalies  with 
a  figure  layout  identical  to  that  used  in  the  corresponding  reflectrix 
displays  of  Figs.  27  and  29.  At  first  glance,  it  may  seem  that  these  two 
sets  of  figures  are  identical,  but  closer  examination  will  show  that  there 
are  some  easily  perceptible  differences  between  the  timetrix  and  reflec¬ 
trix  presentations.  This  is  particularly  true  in  the  details  of  the  dis¬ 
turbance  of  the  ray  structure  by  the  anomalies.  The  figures  are  shown 
here  primarily  as  an  aid  to  specialists  who  may  otherwise  have  wondered 
how  much  error  is  introduced  in  the  reflectrix  construction  as  a  means  of 


Fig.  30.  THE  TIMETRIX  STRUCTURE  IN  THE  AMBIENT  IONOSPHERIC  MODEL. 
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a.  With  the  local  anomaly 


b.  With  the  outward  tilt 


c.  With  the  inward  tilt 


Fig.  31.  TIMETRIX  STRUCTURES  IN  THE  LOCAL  AND  TILTED  ANOMALIES. 

measuring  time  delay  through  the  Dreit  and  Tuve  theorem.  It  appears  rea¬ 
sonable  to  conclude  that  the  reflectrix  will  serve  most  practical  purposes 
quite  well,  but  that,  for  some  specialized  applications,  the  timetrix  (or 
one  of  several  other  possible  modifications)  might  be  superior 
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n.  With  the  full  cycle  sinusoid  based  above  250  km 


b.  With  the  full  cycle  sinusoid  based  above  200  km 
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c.  Positive  half  cycle  of  sinusoid,  electrons  added  above  210  km 
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d  Negative  half  cycle  of  sinusoid,  electrons  deleted  above  200  km 
Fig.  32  TIMETRIX  STRUCTURES  IN  FHE  CONCENTRIC  ANOMALIES 
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F.  The  "Pa ram  Plot/'  A  Graphical  Construction  Designed  To  Show 
Relationships  among  Rays  ‘ 

It  was  noted  in  the  preceding  section  that  the  reflectrix  (or  time- 
trix)  structures  often  cross  one  another  in  a  way  that  often  makes  it 
impossible  for  the  viewer  to  determine  the  pattern  of  the  lines.  Never¬ 
theless,  the  constructions  do  serve  a  useful  purpose  because  they  show 
the  relationship  of  raypath  range,  group  delay,  takeoff  angle,  frequency, 
and  virtual  height  in  a  concentric  ionosphere.  The  usefulness  of  the  con¬ 
struction  lies  in  the  large  number  of  these  parameters  that  are  included 
on  a  single  graphical  presentation. 

Two  years  ago,  the  author  set  out  with  the  intention  of  finding  the 
optimum  graphical  presentation  of  a  family  of  rays,  which  should  ideally 
show  the  relationship  among  a  variety  of  parameters  by  means  of  curves 
that  do  not  cross  one  another  as  do  the  ref lectrixes.  Among  the  desirable 
ray  parameters  are  the  fo1 lowing: 
radio  frequency 

range,  measure!  along  the  ground 
group  delay 

phase  delay,  or  number  of  wavelengths  along  a  ray 
takeoff  angle  of  the  ray 

arrival  angle  (may  differ  from  takeoff  angle  only  in  a  tilted 
ionosphere) 

apogee  height 

virtual  height,  com: uted  from  takeoff  and  arrival  angle 
virtual  height,  computed  from  group  delay 
absorption  along  the  path 
total  Faraday  turns  along  the  path. 

If  the  geomagnetic  field  is  included  in  the  calculations,  then  of 
course  one  would  like  to  display  all  of  the  above  parameters  for  both  the 
ordinary  and  extraordinary  rays.  Obviously,  the  list  contains  far  too 
many  parameters  to  be  presented  on  a  single  chart.  As  will  be  seen,  the 
readability  of  the  chart  decreases  rapidly  as  the  number  of  parameters 
increases. 
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In  an  effort  to  decrease  the  magnitude  of  the  problem  to  m  iagonble 
dimensions,  interest  was  centered  on  only  a  few  of  the  rav  ->  i  l  aneters : 

(1)  We  neglect  the  magnetic  field  for  this  first  effort. 

(2)  V'e  choose  to  display  range,  group  delay,  and  radio  frequency 
because  these  are  primary  measurable  quantities. 

(3)  We  choose  takeoff  angle  (equivalent  to  arrival  angle  in  a  non- 
tilted  ionosphere)  because  power  density  is  a  function  of  the 
rate  of  change  of  other  parameters  with  respect  to  takeoff 
angle. 

(4)  As  a  fifth  parameter,  we  have  chosen  apogee  height,  that  is, 
the  height  of  the  top  of  any  ray.  This  serves  partly  to  tell 
at  a  glance  the  region  of  the  ionosphere  that  exerts  dominant 
control  over  any  given  ray  since  most  of  the  bending  occurs 
at  apogeq.  Perhaps  most  important,  however,  is  the  fact  that 
it  was  found  by  experimentation  that  the  apogee  could  be  used 
as  one  of  two  axes  in  a  plot  that  had  very  desirable 
characteristics. 

(5)  Since  group  delay  is  almost  proportional  to  range,  lines  of 
constant  group  delay  would  look  very  much  like  lines  of  con¬ 
stant  range.  Furthermore,  comparisons  between  the  time  delays 
of  two  rays  is  generally  carried  out  at  a  fixed  range.  Conse¬ 
quently,  we  have  found  it  desirable  to  subtract  from  each  time 
delay  the  time  it  would  require  for  the  speed  of  light  to  tra¬ 
verse  the  surface  of  the  earth  underneath  the  ray.  We  then  deal 
in  "excess  group  delay,"  Tc ,  which  is  defined  by  the  relation 
Te  =  Tg  “  R/ c  where  Tg  is  group  delay,  R  is  ground  range, 
and  c  is  the  speed  of  light. 


Using  the  five  parameters,  we  plotted  a  sample  set  of  data  in  many 
possible  ways,  using  each  possible  pair  of  parameters  as  abscissa  and 
ordinate.  The  remaining  three  parameters  were  used  to  define  families 
of  curves  plotted  in  the  field.  Thus  there  were  ten  possible  plots  con¬ 
sidered;  of  these,  the  best  was  the  plot  in  which  ground  range  served  as 
the  abscissa  and  apogee  height  served  as  the  ordinate.  On  this  display, 
if  the  ionosphere  is  nontilted,  none  of  the  curves  within  a  family  ever 
cross  one  another;  therefore,  it  is  not  possible  to  make  an  ambiguous 
graph.  Furthermore,  the  three  families  of  parametric  curves,  defined 
by  takeoff  angle,  radio  frequency,  and  excess  time  delay,  all  tend  to 
cross  one  another  at  roughly  60  angles,  which  is  optimum  for  least  con- 
•usion.  finally,  the  overall  shape  of  such  a  plot  is  roughly  such  that 
it  fills  most  of  a  rectangular  area  and  thus  the  information  can  be  com- 
onetly  displayed  on  a  page.  This  latter  criterion  is  important  and 
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failure  to  meet  it  was  one  of  the  most  serious  shortcomings  of  an  other¬ 
wise  desirable  display  in  which  range  and  takeoff  angle  served  as  the 
two  axes . 

Figure  33  shows  a  param  plot,  which  has  been  drawn  with  rays  calcu¬ 
lated  through  the  ambient  ionospheric  model.  A  single-color  param  plot 
is  extremely  difficult  to  read,  partly  because  the  computer-drawn  curves 
have  a  small  scale  jitter  that  serves  to  break  up  the  flow  of  the  line 
pattern.  This  makes  it  difficult  for  the  eye  to  follow  the  progress  of 
a  family  of  lines  in  the  presence  of  other  similar  families.  It  was 
possible  to  read  the  single-color  param  plot  in  the  ambient  ionosphere, 
but  even  the  addition  of  a  simple  E  layer  to  the  model  introduced  enough 
additional  complexity  to  make  the  single-color  version  virtually  useless. 
This  is  probably  the  main  disadvantage  of  the  param  plot  representation, 
since  multicolor  reproduction  like  that  of  Fig.  33  is  expensive  in  these 
reports  and  virtually  impossible  to  incorporate  in  other  media  For  ex¬ 
ample,  technical  journals  rarely  reproduce  color,  and  no  office  copying 
machines  will  preserve  the  color  of  a  picture. 

To  overcome  the  need  for  color,  the  remaining  param  plots  in  this 
report  will  be  presented  with  only  the  curves  that  are  parametric  in  fre¬ 
quency.  This  will  serve  our  purpose  here,  but  does  not  solve  the  basic 
problem.  It  would  be  very  useful  if  a  practical  parameter  plotting  for¬ 
mat  could  be  devised,  for  then  the  output  'f  a  raytracing  program  could 
be  much  more  readily  preserved  and  communicated  to  other  workers.  Figure 
33  is  included  here  primarily  to  show  what  we  are  striving  for,  in  the 
hope  that  other  workers  may  become  interested  in  this  display  problem  and 
attempt  to  improve  upon  it. 

In  a  sense,  a  param  plot  is  a  map  of  the  characteristics  of  a  par¬ 
ticular  ionospheric  profile  as  it  affects  radio  waves.  So  long  as  all 
wave  sensors  are  on  the  earth's  surface,  the  param  plot  presents  most  of 
the  useful  Information  that  can  be  obtained.  For  example,  the  radio  sig¬ 
nal  strength  is  a  direct  and  fairly  simple  function  of  the  derivative 
dp^R,  which  can  be  obtained  directly  from  the  param  plot.  Although  not 
presented  here,  phase  path  difference  information  is  inherently  Included 
in  the  form  of  excess  group  delay.  Thus,  one  could  postulate  a  change  in 
critical  frequency  and  thereby  compute  the  induced  frequency  shift  in  some 
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Fig.  33.  THREE  COLOR  PARAM  PLOT  IN  THE  AMBIENT  IONOSPHERIC  MODEL. 


hypothetical  carrier  signal.  Possibly  the  only  serious  deficiency  is  the 
lack  of  absorption  or  Faraday  rotation  information  on  the  chart,  but  i 
is  clearly  impractical  to  try  to  add  other  families  of  curves  to  Fig.  33. 
One  might  consider  supplying  additional  transparent  overlays  providing 
additional  information  on  the  same  abscissa  and  ordinate  axes;  this  should 
work  well.  However,  from  a  practical  standpoint,  ic  is  even  more  objec¬ 
tionable  than  the  use  of  color  to  require  the  availability  of  transparent 
overlays. 

Figure  34a  shows  the  frequency  curves  of  a  param  plot  generated  from 
the  rays  computed  in  the  outward  and  inward  tilt  models.  These  curves 
should  be  compared  with  the  red  curves  on  Fig.  33;  it  can  be  seen  that 
there  is  indeed  a  disturbance  with  a  tilt  corresponding  to  the  anomaly 
tilt.  It  is  the  author's  impression  that  most  of  the  small  structural 
detail  should  be  ignored  because  it  is  orobably  inserted  by  the  imperfect 
ray  calculation  method  that  was  used.  However,  a  very  important  point 
can  be  made  on  this  drawing  which  has  not  been  evident  in  any  of  the  pre¬ 
vious  data  displays,  and  therein  lies  the  value  of  the  param  plot. 

The  tilted  vertical  bars  in  the  drawings  enclose  the  region  contain¬ 
ing  most  of  the  rays  that  have  their  apogees  within  the  disturbances.  It 
is  only  possible  to  draw  such  a  region  on  the  plot  because  of  the  choice 
of  axis  labels,  range,  and  height.  If  the  ionosphere  had  contained  no 
tilts  whatsoever,  then  each  ray  apogee  would  be  located  at  exactly  half 
of  the  total  range  of  the  ray.  With  this  idea  as  an  approximation,  the 
tilted  vertical  bars  are  drawn  by  plotting  the  height  vs  range  of  the 
inner  and  outer  surfaces  of  the  anomalies,  allotting  a  factor  of  two  to 
allow  for  the  fact  that  the  abscissa  represents  a  total  horizontal  range 
rather  than  apogee  range.  Thus,  it  is  not  strictly  true  that  the  bounda¬ 
ries  include  all  those  rays  and  only  those  rays  whose  apogees  lie  in  the 
anomalies  because  the  rays  are  not  symmetric  in  this  ionosphere.  (This 
follows  because  the  ionosphere  is  not  concentric  with  the  earth.)  Never¬ 
theless,  it  can  be  seen  that  the  primary  effect  of  the  anomalies  is  exert¬ 
ed  only  on  those  rays  that  fall  within  the  boundaries  shown  on  Fig.  34a 
and  b. 

It  can  be  seen  why  this  should  happen  if  one  applies  the  logic  that 
was  described  in  Chapter  V  during  the  discussion  of  the  ray  drawings 
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through  the  outward  and  inward  tilts.  The  effect  is  simpie  when  one 
examines  the  disturbance  of  the  lower  rays  only.  Those  rays  that  strike 
the  anomaly  on  their  upward  transit  acquire  a  shortened  range  and  those 
rays  that  strike  it  on  their  downward  passage  acquire  a  lengthened  range 
so  that  the  meeting  of  the  two  families  takes  place  near  the  rays  with 
apogees  in  the  disturbance.  This  simplicity  of  logic  can  be  seen  in  the 
param  plots  if  one  examines  the  lower  rays  from  20  MHz  downward.  Each 
ray  family  appears  comparatively  unaffected  except  at  the  location  where 
the  apogees  pass  through  the  anomaly,  where  there  is  a  noticeable  dis¬ 
continuity.  However,  notice  on  Fig.  34a  that  the  24-  and  26-MHz.  effect 
is  comparatively  complicated  since  it  occurs  near  the  skip  distance  where 
the  anomaly  can  change  upper  rays  to  lower  and  vice  versa.  The  inward 
tilt  causes  this  effect  at  22  and  24  MHz. 

Now  it  can  clearly  be  seen  why  the  height-range  behavior  of  the  anom¬ 
aly  affects  the  frequency-time  behavior  of  the  ground  backscatter.  The 
excess  time  delay  changes  from  about  ,14  msec  at  6  MHz  to  .45  msec  at  24 
MHz  as  we  follow  the  progress  of  the  anomnly  effect  upward  through  the 
rays  on  the  param  plots.  (These  curves  have  not  been  included  here,  but 
they  were  computed. )  For  the  purpose  of  this  discussion,  we  can  assume 
that  the  rr.nge  axis  also  represents  the  time  axis  since  this  change  in 
excess  deity  is  comparatively  minor.  The  relation  between  height  and  fre¬ 
quency  is  almost  linear  as  can  be  seen  by  examining  the  spacing  between  the 
various  frequency  curves  at  their  intersection  with  the  tilted  anomaly. 
Thus,  altitude  is  equivalent  to  frequency  during  the  passage  of  the  radio 
waves.  Similarly,  range  is  converted  into  time  delay  because  the  total 
distance  exercises  dominant  control  over  the  time  delay.  The  ionosphere 
adds  some  effect,  tending  to  increase  time  delay  at  higher  frequencies, 
but  it  is  a  minor  addition  to  the  effect  of  the  tilt.  This  is  important 
in  understanding  the  reason  for  the  upward  tilt  bias  in  the  backscatter 
streaks,  however.  Refer  to  Fig.  8  and  notice  that  the  local  anomaly  pro¬ 
duces  a  streak  that  is  tilted  upward  about  1°.  The  outward  tilt  adds  a 
positive  2"  increment  and  a  negative  tilt  subtracts  2“  from  this  1°  pre¬ 
vailing  tilt.  The  prevailing  tilt  is  attributable  to  the  effect  of  the 
ionosphere.  The  additional  2°,  plus  or  minus,  is  attributable  to  the 
tilt  of  the  anomaly 
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Thus,  the  param  plots  have  permitted  us  to  deduce  the  reason  for  the 
1°  bias  in  the  tilts  on  Ftp.  8. 

Figure  35  shows  param  plots  for  the  four  concentric  anomalies.  Again, 
it  is  reasonably  valid  to  assume  that  the  range  axis  provides  a  measure  of 
group  delay  since  the  excess  delay  curves  add  a  comparatively  minor  pertur¬ 
bation  to  numbers  computed  by  simply  converting  range  values  to  time  val¬ 
ues.  Figure  35a  shows  the  frequency  curves  for  the  ravs  used  to  generate 
the  backscatter  of  Mg.  8d.  The  concentric  anomaly  was  located  high 
enough  in  the  ionosphere  that  there  were  no  propagating  rays  with  apogees 
in  the  ambient  model  above  the  anomaly.  This  explains  why  the  backscatter 
from  the  high  concentric  anomaly  was  more  simple  than  the  backscatter  from 
the  lower  concentric  anomaly.  Figure  35b  shows  the  lower  anomaly  nnd  now 
we  see  that  there  are  a  large  number  of  rays  with  apogees  above  the  anom¬ 
aly.  Most  of  the  single  frequency  curves  have  three  distinct  minimum  range 
points.  (They  similarly  have  three  minimum  time  points.)  Consequently, 
i here  should  be  three  distinct  enhancements  in  any  single  frequency  back- 
scatter  calculation;  indeed,  this  was  the  case  in  Fig.  20. 

On  Fig.  35b  notice  that  the  upper  two  noses  on  any  single  curve  are 
never  far  apart  in  range.  This  is  why  the  added  enhancement  in  the  back¬ 
scatter  never  departed  very  much  from  the  normal  leading  edge  of  the  back¬ 
scatter  on  Fig.  20a.  However,  notice  that  the  comparison  of  Fig.  35n  and 
b  shows  a  systematic  difference  between  the  structure  of  the  rays  nnd  the 
two  ionospheres.  The  nmbient  leading  edge  of  the  backscatter  is  preserved 
on  Fig.  35a  in  the  form  of  rays  that  reach  their  minimum  range  at  an  alti¬ 
tude  of  roughly  230  km,  well  below  the  anomaly.  However,  on  Fig.  35b,  the 
reproduction  of  the  ambient  leading  edge  of  the  backscatter  is  attributable 
to  rays  that  have  their  apogees  at  about  250  to  270  km  well  above  the  con¬ 
centric  anomaly,  Presumably,  an  anomaly  located  at  some  intermediate  alti¬ 
tude  would  cause  a  different  effect.  Thus  it  appears  that  a  traveling 
anomaly  moving  vertically  in  the  ionosphere  should  cause  widespread  dis¬ 
turb  mces  throughout  the  body  of  any  sweep  frequency  backscatter  echo 
structure.  It  should  be  possible  to  gain  a  full  understanding  of  the  phe¬ 
nomena  by  using  ray  plots,  param  plots,  and  synthetic  backscatter  plots 
in  conjunction  with  each  other.  We  have  not  carried  this  study  further, 
but  it  does  appear  to  be  an  area  in  which  more  research  is  desirable. 
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a.  With  the  sinusoid  above  250  km 


Figure  35c  and  d  show  the  effect  of  propagation  in  the  partial  con¬ 
centric  anomalies  which  are  generated  by  separating  the  positive  and  nega¬ 
tive  half  cycles  of  the  sinusoid  disturbance  which  was  used  to  generate 
ic nosphere  I  ID  319.  The  reason  for  the  structure  of  the  backscatter  can 
now  be  seen  by  applying  logic  similar  to  that  discussed  above.  It  does 
appear  on  the  basis  of  this  calculation  that  it  would  be  possible  to  make 
a  simple  electron  addition  of  such  a  nature  that  the  third  backscatter  en¬ 
hancement  is  suppressed,  but  it  would  require  careful  model  structuring. 

Perhaps  the  lesson  to  be  learned  is  that  the  appearance  of  two  new  back¬ 
scatter  enhancements,  before  and  after  the  ambient  leading  edge,  are  in¬ 
dications  of  a  single  phenomenon  similar  to  that  shown  in  Fig.  35b,  c, 
or  d. 
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VII  CONCLUSIONS 


In  displays  of  sweep  frequency  backscatter  sounding  data  there  are 
often  streaks  in  the  data  which  are  nearly  parallel  to  the  frequency  axis. 
Here  it  has  been  shown  that  these  streaks  may  be  attributed  to  the  pres¬ 
ence  of  localized  irregularities  in  the  ionosphere  approximately  midway 
between  the  sounder  and  the  scattering  ground.  When  striations  nre  tilt¬ 
ed  (that  is,  when  their  time  delay  varies  with  frequency),  the  tilt  may 
be  attributed  to  a  geometric  tilt  in  the  parent  anomaly.  This  cause- 
effect  relationship  has  been  partially  verified  by  an  observation  that 
the  direction  of  tilt  and  the  direction  of  motion  of  individual  striations 
agree  with  the  predictions  of  atmospheric  gravity  wave  theory,  at  least 
for  the  six  examples  available.  It  was  also  found  that  the  frequency  in¬ 
terval  spanned  by  a  backscatter  striation  should  indicate  the  height  in¬ 
terval  occupied  by  the  anomaly. 

Finally,  it  has  been  shown  that  horizontally  disposed  anomalies 
(  layerlets ")  in  an  otherwise  smooth  ionospheric  layer  give  rise  to  new 
leading  edges  on  the  backscatter  records.  Other  mechnnisms  enn  enuse 
similar  multiple  lending  edges;  discrimination  among  the  possible  causes 
may  be  feasible  but  has  not  been  demonstrated. 
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APPENDIX 


FIG.  3.  BACKSCATTER  AT  2  MINUTE  INTERVALS,  ARRANGED  TO  SHOW  MOTION  BY 
PAGE  FLIPPING. 
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Hfy  use  of  computer  simulation,  it  is  possible  to  predict  the  effect  of  traveling 
ionospheric  Irregularities  upon  the  data  obtained  by  HF  sweep-frequency  ground- 
backscatter  sounders.  The  simulation  process  makes  use  of  an  ionospheric  model  in 
which  radio  rays  are  traced  to  determine  the  time  distribution  of  energy  along  the 
distant  ground  after  skywave  travel.  Subsequent  computation  yields  a  synthetic  record  < 
of  the  echo  amplitude  as  a  function  of  radio  frequency  and  time  delay,  much  like  the 
experimental  data.  The  degree  of  realism  achieved  is  limited  primarily  by  the  need  for 
economy;  the  process  used  here  is  a  simplified  version  of  one  previously  described 
(Rei.  1). 

It  is  shown  that  a  localized  ionospheric  irregularity  causes  a  discontinuity  in 
backscatter  echo  amplitude  with  a  limited  extent  in  frequency  and  with  a  delay  that  is 
comparatively  independent  of  frequency.  This  discontinuity  has  the  appearance  of  a 
short  streak  on  the  synthetic  records;  such  a  feature  is  common  in  experimental  data. 
When  an  ionospheric  irregularity  exists  from  the  ground  to  an  unlimited  altitude  and  is 
formed  with  its  long  axis  tilted  about  45°  away  from  the  observer,  a  backscatter  streak 
appears  which  is  spread  over  the  entire  frequency  range  of  the  background  echo.  Fur¬ 
thermore,  this  streak  is  tilted  on  the  record,  i.e.,  its  time  delay  increases  with  an 
increase  in  frequency.  When  a  simile  anomaly  is  programmed  so  that  it  tilts  inward 
toward  the  observer,  then  the  resulting  streak  has  a  delay  that  decreases  with  an  in¬ 
crease  in  frequency.  A  relation  between  the  tilt  of  the  anomaly  and  the  tilt  of  the 
backscatter  streak  is  thus  demonstrated.  The  relation  is  in  accord  with  a  limited  ob¬ 
servation  that  tilted  streaks  in  experimental  backscatter  move  only  in  the  predicted 
direction  of  anomaly  motion  as  deduced  from  the  measured  direction  of  anomaly  tilt. 
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Abstract  (contd) 

Anomalies  with  constant  altitude  are  shown  to 
produce  new,  short-range  leading  edges  on  the  back- 
scatter.  Experimental  data  having  this  form  can 
thus  be  interpreted  as  an  indication  of  the  presence 
of  small  "layer lets"  within  the  ionosphere,  but  it 
is  noted  that  there  are  at  least  three  other  pos¬ 
sible  causes  for  this  particular  form  of  backscntler 
echo. 


